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ABSTRACT 


A  single  dimensional  sinusoidal  parameter  perturbation  adaptive 
control  system  Is  studied.  The  operation  and  signals  at  every  point  of 
the  system  are  analyzed  theoretically.  Then,  the  system  Is  simulated  on 
analog  computer  starting  with  the  simplest  circuit  and  adding  more  com¬ 
ponents  until  a  complete  system  Is  formed.  The  results  of  simulation 
study  are  compared  with  the  theoretical  analysis.  Both  agree  In  most 
cases.  Finally,  an  analytical  design  of  the  system  Is  attempted  based 
on  the  simulation  study  and  using  a  linearizing  technique.  It  was  simu¬ 
lated  on  analog  computer  and  the  performance  was  compared  with  the  one 
predicted  analytically.  It  does  not  follow  exactly  as  predicted  but  the 
deviation  Is  not  so  great.  This  method  Is  useful  In  obtaining  preliminary 
design  which  will  be  refined  by  experimental  verifications. 

The  writer  wishes  to  express  his  appreciation  for  the  assistance  and 
encouragement  given  him  by  Professor  R.  C.  Dorf  of  the  U.  S.  Naval  Post¬ 
graduate  School  in  this  Investigation. 
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1.  Introduction 


The  adaptive  control  system  has  received  much  attention  recently. 

The  principle  of  adaptive  control  is  to  maintain  a  system  always  in  an 
optimum  operating  condition  whatever  the  affecting  situation  may  be.  There 
are  two  kinds  of  adaptive  control  systems;  the  signal  adaptive  system  and 
the  system  adaptive  system.  The  signal  adaptive  system  is  a  system  which 
can  analyze  incoming  signals  and  automatically  adjust  the  system  to  give 
an  optimum  output.  The  system  adaptive  system  is  a  system  which  can  mea¬ 
sure  its  own  characteristics  and  then  change  its  structure  or  design  so  as 
to  continuously  maintain  a  desired  performance.  This  paper  deals  with  only 
the  latter.  The  system  adaptive  system  should  have  the  following  two  es¬ 
sential  features  to  achieve  its  goal; 

(1)  The  measurement  of  its  own  characteristics. 

(2)  The  generation  of  proper  corrective  signal  for  the  control¬ 
led  system  according  to  the  measurement.  In  the  measurement  of  the 
characteristics  of  the  system,  the  optimum  operating  condition  or  figure 
of  merit  is  defined.  The  figure  of  merit  or  index  of  performance  (IP.)  is 
an  indication  of  the  characteristic  or  quality  of  a  system.  The  IP  may  be 
a  physically  meaningful  value  such  as  heat  produced  per  gallon  of  fuel  or 
theoretical  quantity  such  as  RMS  error  or  mean  squared  error. 

The  sinusoidal  parameter  perturbation  adaptive  control  system  is  a 
method  of  system  adaptive  control  and  in  this  paper  the  operation  and 
signals  at  every  point  of  this  system,  are  analyzed  theoretically.  Then, 
the  system  is  simulated  on  analog  computer  starting  from  the  simplest  cir¬ 
cuit.  Adding  more  components  one  by  one,  the  effect  of  each  component  is 
observed  until  the  complete  system  is  formed.  The  results  of  these  simula¬ 
tion  study  are  compared  with  the  theoretical  analysis.  Both  agree  in  most 
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cases.  Finally,  a  design  of  the  system  is  attempted  analytically  based  on 
the  simulation  study  and  using  a  linearizing  technique.  It  was  also  simu¬ 
lated  on  the  analog  and  the  digital  computer  and  the  actual  performance  of 
the  system  is  compared  with  the  predicted  performance  from  the  analytical 
design.  The  results  show  that  the  actual  performance  does  not  exactly  agree 
with  the  performance  predicted  analytically  but  it  follows  reasonably  close 
pattern  predicted  analytically.  This  is  expected  because  the  analytical 
method  involves  approximation.  This  method  is  useful  to  obtain  a  prelimin¬ 
ary  design  which  will  be  refined  by  experimental  verification  and  gives  a 
direction  of  improving  the  first  trial  design. 
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2*  The  Sinusoidal  Parameter  Perturbation  Adaptive  Control  System* 

A  single  dimensional  sinusoidal  perturbation  system  is  as  Fig*  1. 

The  model  is  an  optimum  operating  model.  When  the  plant  is  operating 
exactly  as  the  model  there  will  be  no  error  output.  But  when  a  parameter 
in  the  plant  is  changed  to  some  other  from  the  optimum  parameter  value  due 
to  environment  change,  there  will  be  an  error  output.  By  squaring  the  er¬ 
ror  we  get  a  non-negative  function  of  error.  As  an  index  of  performance 
we  may  adopt  average  of  error  squared  as: 


IP 


2 

e  dt 


(2-1) 


This  quantity  is  a  function  of  the  characteristics  of  input  to  the  plant 
and  also  of  the  characteristic  of  plant  parameter  drift  caused  by  environ¬ 
ment  change*  Fig.  2  shows  the  contours  of  IP(X,Y)  for  two  dimensional 
case  where  IP  is  a  function  of  two  parameters  X  and  Y.  The  minimum  IP 
occurs  at  X=XQ  and  Y^Y^.  If  we  fix  the  value  of  Y  at  Y^  and  vary  only  X, 
the  IP  variation  with  X  becomes  as  Fig.  3,  the  single  dimensional  case. 
Generally  the  IP  with  one  parameter  variable  tends  to  be  a  parabolic  func¬ 
tion  of  the  parameter  and  can  be  approximated  near  the  optimum  parameter 
value  as; 

IP(X)  =  ai(X-X0)2  +  a2  (2-2) 

where  X^  is  optimum  parameter  value  and  a^  and  a^ 
are  arbitrary  constants  which  change  the  shape 
of  IP  curve  such  that  it  represents  closely  the 
actual  IP  variation  of  any  plant. 

In  some  plants  a  parameter  has  a  slight  or  no  change  with  respect  to  the 
environment  change  and  such  parameter  is  fixed  at  the  optimum  value  and  the 
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DETECTION  REFERENCE 


(a)  General  set-up 


(b)  Model  Reference  System 


Pig.  1.  Single  Dimensional  Sinusoidal  Parameter  Perturbation  Adaptive 
Control  System. 


Pig.  3.  IP  v.s.  parameter  X 


IP(X,Y) 


Pig.  2.  Contours  of  IP(X,Y)  for  two  dimensional  case. 


IP(X) 
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number  of  parameters  to  be  controlled  is  reduced. 

Referring  to  Fig.  3,  suppose  we  perturb  the  parameter  X  sinusoidally 
by  a  perturbation  signal,  sin  CO  ^ t.  If  the  actual  plant  para¬ 

meter  X  Is  at  the  optimum  parameter  value,  X^,  there  will  be  no  IP  varia¬ 
tion  output  at  frequency  CO  y  If  the  X  value  is  lower  than  X^,  there  will 
be  an  IP  variation  output  at  frequency  CO^  and  the  output  signal  will  be 
180°  out  of  phase  with  respect  to  the  input  perturbation  signal  as  seen  in 
the  Fig.,  on  the  contrary,  if  the  X  is  higher  than  the  X^,  the  output  is  in 
phase  with  the  input  perturbation  signal.  It  is  also  noted  that  the  ampli¬ 
tude  of  the  output  is  proprotional  to  [x-X^j,  therefore  the  output  is  consid¬ 
ered  to  be  a  carrier,  modulated  by  the  variation  of  jx-Xj  through  the  IP 
characteristic.  The  frequency  of  carrier  is  CO ^  and  it  carries  the  informa¬ 
tion  about  the  polarity  of  the  x'xo  by  its  phase  difference  with  the  input 
perturbation  signal.  When  this  output  is  demodulated  by  multipling  with 
the  original  perturbation  signal,  it  will  produce  a  D.C.  signal  with  proper 
sense  to  correct  the  plant  parameter  deviation. 

The  bandpass  filter  eliminates  unwanted  signal  components  contained  in 
the  modulated  signal  and  also  affects  the  response  of  the  adaptive  loop  to 
a  step  change  of  signal. 

The  detector  is  a  multiplier  which  multiplies  the  modulated  signal  with 
original  perturbation  signal  and  produces  positive  or  negative  correction 
signal. 

The  filter  after  the  detector  is  a  smoothing  filter  which  filters  out 
high  frequency  components  in  the  correction  signal  and  acts  as  compensator 
to  stabilize  the  adaptive  loop. 

The  integrator  is  a  D.C.  motor  in  practical  operation  and  it  integrates 
correction  signal  and  drives  the  plant  parameter  control  to  the  optimum  value. 
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If  we  want  to  make  more  than  one  parameter  adaptive,  separate  loops 
for  each  different  parameter  must  be  established  and  use  different  perturba¬ 
tion  frequency  for  each  parameter*  The  same  error  measure  can  be  used  for 
each  loop  if  it  is  suitable,  but  for  each  parameter,  other  elements  must  be 
duplicated. 

Sometimes,  a  plant  has  several  local  minimum  in  its  IP  characteristic. 
In  that  case,  a  special  measure  is  taken  to  find  the  true  minimum.  But 
generally,  all  plants  have  a  single  minimum  around  the  optimum  parameter 
value*  The  sinusoidal  perturbation  method  can  also  be  applicable  to  an 
adaptive  system  which  searches  for  a  maximum  IP  value. 


7 


3.  The  Analysis  of  Signals  in  the  Adaptive  System 
3.1  Block  diagram  of  an  adaptive  system 

To  set  up  a  block  diagram  of  a  single  dimensional  parameter  perturba¬ 
tion  adaptive  control  system,  the  circuit  models  of  each  component  are 
established.  Referring  to  Fig.  1,  the  parameter  perturbation  signal  per¬ 
turbs  a  plant  parameter  and  the  parameter  itself  also  is  changing  its  value 
due  to  environment  change  and  parameter  correction  signal  coming  through 
adaptive  loop.  This  parameter  perturbation  signal  and  parameter  drift 
cause  the  error  output  while  the  plant  and  model  are  operating  with  input, 
R.  It  is  necessary  to  find  out  an  equivalent  transfer  function  block  be¬ 
tween  the  parameter  perturbation  input  and  the  error  output.  The  para¬ 
meter  drift  is  also  an  input  for  the  error  output  and  the  input  is  same  as 
perturbation  input.  We  may  designate  this  equivalent  transfer  function  as 
error  measurement  transfer  function  and  express  as 


K 


m 


G  (s)  =  -  _ 

mK  '  1  +  T  s 


0-1) 


K  and  *7?  are  decided  by  the  signal  transfer  characteristic  between 
perturbation  input  and  error  output.  To  find  out  the  transfer  character¬ 
istic,  phase  lag  and  attenuation  introduced  to  this  error  output  with  res¬ 
pect  to  the  original  perturbation  and  parameter  drift  are  observed.  Al¬ 
though  the  perturbation  and  drift  do  not  come  into  the  plant  through 
ordinary  input  point  but  through  parameter  itself,  these  signals  encounter 
the  same  phase  shift  as  a  signal  of  the  same  frequency  coming  into  the 
plant  through  the  input  of  the  plant,  encounteres.  This  will  be  confirmed 
by  actual  simulation  study  in  Section  4.1.  The  attenuation  factor  of  this 
error  signal  is  not  the  same  as  the  attenuation  factor  of  the  signal  of  same 
frequency  coming  through  plant  input.  The  attenuation  factor  to  the  error 
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signal  depends  on  the  frequency  and  magnitude  of  the  plant  input,  R  and 
the  characteristic  of  the  plant.  It  is  also  governed  by  its  original 
signal  frequency.  If  we  fix  the  amplitude  and  frequency  of  the  input  to 
the  plant,  the  attenuation  factor  is  a  function  of  its  original  perturba¬ 
tion  or  drift  signal  frequency  assuming  the  plant  characteristic  are  con¬ 
stant. 

From  the  above  discussion,  the^  is  chosen  such  that  the  equivalent 

m 

transfer  function  gives  the  same  phase  lag  as  the  plant  introduces  to  the 

same  signal  coming  through  the  plant  input. 

Next,  the  K  is  decided  by  actual  measurement  of  attenuation, 
m 

In  measuring  the  attenuation,  the  input  to  the  plant  is  kept  constant  and 
perturb  or  drift  the  parameter  by  some  signal  and  measure  the  error  out¬ 
put,  which  gives  the  attenuation  factor  between  parameter  drift  input  and 
the  error  output.  Let  the  measured  attenuation  be  A  .  Then 


A 

m 


) 


(3-2) 


where  (4)^  is  a  parameter  drift  frequency  used 

for  measuring  A  . 

m 

From  equation  (3-2)  K  can  be  calculated. 

m 

2 

The  IP  curve  is  represented  by  IP(X)  =  a^(X-X^)  +  as  explained  in 

Section  2,  A  multiplier  is  used  as  squarer  and  a^  is  the  coefficient  of 
the  squarer.  After  the  squarer,  is  added.  A  proper  choice  of  a^  and 
gives  an  approximate  representation  of  any  single  dimensional  IP  curve. 

The  bandpass  filter  has  the  following  form, 

CO2  s 

G.(s)  =  T - - -  2  (3-3) 

1  S  +  2/^S  +  CO 

The  detector  M,  is  a  multiplier. 
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The  smoothing  filter  ta^es  various  forms  depending  on  the  re¬ 

quirement  of  compensation  of  the  adaptive  loop  but  generally  it  is  a  phase 
lag  filter. 

The  reference  signal  filter  is  an  arbitrary  filter  chosen  such  that 
phase  shift  at  frequency  due  to  this  filter  is  equal  to  the  phase 
shift  introduced  to  the  perturbation  signal  between  the  perturbation  in¬ 
put  and  detector  through  the  adaptive  loop. 

The  integrator  is  usually  a  D.  C.  motor  and  its  transfer  function  is 
designated  as  G^(s), 

The  complete  block  diagram  of  a  single  dimensional  parameter  perturba¬ 
tion  adaptive  control  system  becomes  as  Fig.  4. 


10 


11 


Fig.  4  Block  diagram  of  a  single  dimensional  sinusoidal 
parameter  perturbation  adaptive  control  system. 


3.2  Modulation 


\ 


The  input  signals  through  the  parameter  are  the  perturbation  signal 
and  parameter  drift  signal  due  to  the  variation  of  the  parameter  value. 
The  instant  error  deviation  of  parameter,  Y^(t)  can  be  expressed  as 

YjCt)  =  X(t )  +  Cp  -  XQ  (3-4) 

where  X(t)  =  instantaneous  parameter  drift 
X0  =  optimum  parameter  value 

=  parameter  perturbation  signal 

The  analysis  of  signals  will  be  made  for  two  cases  of  parameter  change; 
sinusoidal  drift  of  parameter  and  step  drift  of  parameter. 

(a)  Sinusoidal  drift  of  parameter 

The  parameter  drifts  sinusoidally  around  the  optimum  parameter  value 
with  very  low  frequency,  =  C^sin^t  is  used  as  perturbation  signal 

throughout  the  analysis.  The  instantaneous  parameter  drift  is 

X(t)  =  X^  +  X-  sin  CO  t  (3-5) 

0  1m 

00  -  parameter  drift  frequency  and  is  much  lower 

m 

frequency  than  perturbation  frequency, 

Then  Y^(t)  becomes  as  below  using  equation  (3-4), 

Y-(t)  =  C-sin  66>  t  +  X,  sinft^  t  (3-6) 

1  111m 

Referring  to  the  Fig.  4,  Y^(t)  goes  through  G^(s)  and  comes  out  as 
Y2(t), 

Y„(t)  =  e(t)  =  A  C .  s  in  (  CO  t  -  0  )  +  A  X,sin(^  t  -  0  )  (3-7) 

Z  pi  lpml  mm 

where  A  ,  A  are  associated  attenuation  factors,  and 
p  m 

9  ,  0  are  associated  phase  shifts. 

p  m 
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Signals  at  each  point  in  the  loop  are  labeled  in  the  Fig.  A. 
Following  the  loop, 


Y3(t)  =  a 


1  [V'>] 


Y4(t)  =  a2  +  Y3(t) 


Calculating  above  equations,  Y^(t)  becomes 


Vc> 


a2c2, 

a2  +  ai^  2 


2  2  2  2 
A  X.  a  A  C.  0 

m  1 _ *  1  p  1  cos2 

*2  ;  2 


2  2 

aiAm"l  cos2  (CO  t  -  9  ) 
- 2 -  m  m 


(3-8) 

(3-9) 


<«Y  -  V 

(3-10) 


+  2a.A  A  C-X,sin(&)  t  -  9  )sin(&/,t  -  9  ) 

•  X  tn  p  1  1  m  m  1  p 

The  last  term  of  Y,(t),  2a,A  A  C,X,sin(^  t  -  9  )sin(tc>,t  -  9  )  is  the 

4  lmpll  m  nr  v  1  p 

modulated  signal  term  where  the  perturbation  signal  is  amplitude  modulated 
by  parameter  drift  signal.  This  term  carries  not  only  the  information  of 
magnitude  of  drift  signal  but  also  the  information  about  the  actual  operat¬ 
ing  point  of  parameter,  in  other  words,  the  actual  parameter  is  higher  or 
lower  than  the  optimum  parameter  value  by  the  following  form  in  an  actual 
signal . 

When  X(t)  -  XA  )  0  or  X1sin(&?  t  -  9  )  >  0,  this  term  is 
0  1  mm 

2a  A  A  C,  I  X.  sin(  t  -  9  )  I  sin(aAt  -  9  )  (3-11) 

Impl,I  m  mi  i  p 

When  X(t)  -  0  or  Xtsin(&>  t  -  9  )/  0,  this  term  is 

7  0  1  v  m  m  ^  * 

2a  A  A  C.  I  X-sin(&>  t  -  9  )  I  sin(cO,t  -  9  -  rt)  (3-12) 

lmpl|  1  m  nr  I  v  1  p  7  7 

This  is  explained  in  Section  2  and  shown  in  Fig.  3. 
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(b)  Step  drift  of  parameter 

The  parameter  value  drifts  in  step  from  optimum  parameter  value  to 
some  value  X,  then  by  equation  (3-4),  the  instantaneous  deviation  of  para¬ 
meter  becomes 

Y.(t)  *  X  +  C  -  Xn  »  (X  -  O  +  C.sin&At  (3-13) 

1  p  U  U  1  1 

Going  through  the  equivalent  transfer  block,  G  (s)  and  comes  out  as 

m 

Y2(t)  =  e (t)  =  Ad(X  -  XQ)  +  ApCjSinC^t  -  9  )  (3-14) 


where  A.,  A  are  associated  attenuation  factors  and 
d*  P 

0  is  the  phase  shift. 

P 

Calculating  the  same  way  as  in  case  (a)  using  equations  (3-8)  and  (3-9), 


Y^(t)  comes  out  as 


Y4(t)  *  a2 


+  a 


A2(X- 


xo)2  + 


aV 

p  1 


■) 


.  2  2 
a  A  C 

1  ?  i 


COStcl^t 


+  2a  A  A  C.(X  -  X^)sin(&J.t  -  0  )  (3-15) 

1  d  p  1  0  1  p 

In  Y,  (t)  the  last  term,  2a. A  A  C.  (X  -  Xrx)sin(^;it  -  0  )is  an  amplitude  modulated 
4  7  ’  1  d  p  1  01  p7  r 

signal  term  where  the  carrier  is  the  perturbation  signal,  sin  (ftAt  -  0  ). 

1  P 

This  term  carries  the  information  about  the  amount  of  parameter  drift  (X-X^) 
and  furthermore  the  direction  of  parameter  drift  from  the  optimum  parameter 
value.  This  term  takes  following  form  for  actual  signals: 

When  X  -  X^  0  or  the  actual  parameter  is  higher  than  the  optimum 
parameter  value,  this  term  becomes 

2a. A.A  C.  I  X  -  X„  I  sin(&>.t  -  0  )  (3-16) 

ldpll  0  '  1  p7  v  7 

When  X  *  X^  0  or  the  actual  parameter  is  lower  than  this  becomes 
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2aiVpci  I  x  -  x0  I  -  ep  -  t) 


as  explained  in  Section  2  and  shown  in  Fig.  3. 


(3-17) 
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3o  3  Filtering 

The  bandpass  filter  G^(s)  attenuates  considerably  all  unnecessary 
signals  in  Y^(t)  and  only  passes  the  components  that  contain  information 
about  the  parameter  drift*  without  much  attenuation. 

(a)  Taking  the  case  of  sinusoidal  parameter  drift*  the  Y^(t)*  equation 
(3-10)  in  Section  3.2(a),  is  rewritten  as 


Y,  (t)  =  D  -  J),cos2(CV.t  -  6  )  as  Dncos2 (iO  t  -  0  ) 
4  0  1  1  p  2  m  nr 


where 


Do  = 


+  D0sin (to  t  -  0  )sin(tJ1t  -  G  ) 
3  mm  1  p 


1  ?  o  o  o 

a2  +  — <Vl  +  W 


D,  = 


,2  2 
a  A  C 

1  p 1 


(3-18) 


D„  = 


2  2 
a  A  X, 
lml 


D-  =  2a  A  A  C.X, 

3  1  m  p  1  1 

Expanding  the  last  term  of  equation  (3-18), 


Y4(t)  =  DQ  -  D1cos2(to1t  -  0  )  -  D2cos2(^mt  -  GJ 


3  cos  r 

2  L 

6») 


( a).  -  io  )  t  -  9  +  e 

1  m  pm 


]  --^  'osh 


+  U)  )t  -  G 

m  p 


(3— 19) 


The.  Y^(t)  contains  the  following  frequency  signal  componentss 

D.C.  s  2U>  ,  2  CO  ,  CV  -O)  ,  CO  +  lO  . 

1  m’  11  m  l  m 
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The  necessary  intelligence  is  contained  in  CO.  -  U)  and  CO .  +  U)  com- 

1  m  1  in 

ponents,  Therefore  the  bandpass  filter  is  designed  to  pass  these  two 
components  without  attenuation  and  suppress  other  components.  Thus  the 
center  frequency  of  the  bandpass  filter  should  be  CO '  and  this  filter  does 
not  introduce  any  phase  shift  to  the  signal  of  cOy  Assuming  that  the 
phase  angle  characteristic  is  linear  near  the  center  frequency*  which  is 
generally  true  for  an  actual  filter,  the  signal  of  (  CO ^  gets  +  A  0 

phase  shift  and  the  signal  of  (  CO  ^  gets  -  A  0  phase  shift  with 

respect  to  the  signal  of  Cl) j.  Then  Y^(t)  passing  the  filter,  becomes  as 


Y5(t)  *  Vo  -  VlC052  <“'lt  •  9p  -  V  •  A2D2  cos2(wmt  -  8b  -  92) 


A.D, 

3  3  cos 


f  (£0  -  (o  )t  -  6  +  0  +  A  0 1 

L  1  m  pm  J 


A.D. 

4  3  cos 


[<"1  +  a> 


)t  -  0  -  0  -  A  0 

m  pm 


(3-20) 


where  A,.  A„,  A0,  and  A.  are  attenuation  factors  associated  eith  each 

0  1  2  3  4 

signal  and 

0^,  and  A  0  are  additional  phase  shifts  introduced  to  each  signal. 
Usually  in  a  band  pass  filter  A^  and  A^  are  the  same  attenuation  as  the 

center  frequency  signal  suffers  because  as  assumed  in  Section  3,2, 

Aq,  A^  and  A2  should  be  much  smaller  than  A^  for  our  purpose.  Combining 
the  last  two  terms  of  equation  (3-20),  ^(t)  becomes 


Ys(t)  =  AqD0  -  A1D1cos2(fc>1t  -  0p 


0.)  -  A0D0cos2(  U)  t  -  0  -  0O) 

X  z  z  m  m  z 


+  A_D_sin(  CO.  t  -  0  )sin(tu  t  -  0  -  A  0) 

3  3  1  p  mm  ' 


(3-21) 
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If  the  modulated  signal  term  in  Y^(t),  equation  (3-18)  is  compared  with 
the  modulated  signal  in  Y^(t),  equation  (3-21),  it  is  noted  that  the  band¬ 
pass  filter  introduced  a  phase  lag  of  A  9  only  to  the  modulating  signal 
component,  (sin  U) ^ t  -  0^)  as  far  as  the  modulating  signal  term  is  concerned. 

(b)  In  step  parameter  drift  case,  Y^(t),  equation  (3-15)  is  rewritten  as 


V° 


D5  -  DgCos2  CO  1 1  +  D7(X  -  XQ)sin(ic»1t 


V 


(3-22) 


where  Dr  = 


a2  +  a 


i  [  Ad(x  -  V2  *  Hr~) 


D,  = 


.22 
3. A  c 

i  p  i 


D_  =  2a  A  A  C, 

7  1  d  p  1 

Passing  the  bandpass  filter,  Y^(t)  becomes 

Y5(t)  =  AqD5  -  AjDgCos (2  CO ^  -  02)  +  A^^X  -  X^sin^jt  -  Gp)  (3-23) 
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3.4  Demodulation 


The  modulated  signal  Y,.(t)  reaches  the  detector  Mg,  where  Y^(t)  is 
multiplied  by  the  detection  reference  signal  and  demodulation  is  accomplish¬ 
ed.  The  detection  reference  signal  should  be  in  the  same  phase  as  the  car¬ 
rier  signal  to  get  maximum  output  at  the  demodulation.  For  this  purpose, 
detection  reference  signal  phase  lag  filter  G^(s)  is  used  to  give  the  refer¬ 
ence  signal  an  exactly  same  phase  lag  as  the  phase  lag  introduced  to  the 
carrier  signal  between  the  perturbation  input  at  the  plant  and  the  detector; 
in  our  case  which  is  0^.  Thus  the  detection  reference  signal  has  the  follow¬ 
ing  form, 

Y6(t)  -  BC^sinC  COl t  -  0p)  (3-24) 

where  B  =  attenuation  factor  of  G^(s) 

0^  =  phase  lag  introduced  by  G^(t) 

(a)  For  the  sinusoidal  parameter  drift  case,  Y,.(t)  equation  (3-21)  is 
multiplied  by  Y^(t)  equation  (3-24)  and  we  get 


Yy(t)  =  A0D0BClSin(  (V  f  -  0  )  -  B  C1cos2(^1t  -  0p  -  0^ 

sin(  Co  t  -  0  )  -  A  D0B  C.cos2(^  t  -  0  -  0o)sin(uAt  -  0  ) 

1  p  z  z  l  m  m  z  lp 

2 

+  AnD_B  C.  s  in  (UJ-t  -  0  )sin(^  t  -  0  -  A0) 

331  v  1  p  m  m  ' 


(3-25) 


In  Y^(t)  only  a  D.C.  component  is  integrated  by  the  integrator  following  and 
develops  a  parameter  correction  signal.  Only  the  last  term  of  Y7(t)  contains 
a  D.  C.  component  which  is  varying  slowing  in  accordance  with  the  modulat¬ 
ing  signal,  which  is  the  parameter  drift.  All  other  terms  do  not  contribute 
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to  the  parameter  correction  because  they  are  comparatively  high  frequency 
sinusoids  and  furthermore,  A^,  and  A^  are  very  small  due  to  the  band¬ 
pass  characteristic  of  G^(s)  and  they  will  be  further  attenuated  by  the 

smoothing  filter  following*  If  only  the  last  term  of  Y-(t)  is  expanded 

7 

it  becomes 

A  T)  u  c  A  D..B  C- 

3  3  Isin  {to  t  -  9  -49)-  - - - sin(  to  t  -  9  -A0) 

-  mm  2  mm 


cosHcQ^t  -  0p)  (3-26) 

The  first  term  of  (3-26)  is  considered  as  a  slowly  varying  D.  C.  which 
corrects  the  parameter  drift.  Thus, 

A  D  B  C 

Y  (t)  =  - COt  -  9  -  4  9)  (3-27) 

/  a  *  c .  Z  m  m 

Putting  back  the  original  value  for  and  taking  into  account  the  phase 
shift  between  Y^(t)  and  Y^(t)# 


Yn(t)  s=  a, A  A  A.B  C,X,sin(  CO  t  -  9  -  A  9  -  9.) 

9  1  m  p  4  11  m  m  3 


(3-28) 


where  0^  is  the  phase  shift  in  ^(s)  anc*  G^(s)* 

A.  is  the  attenuation  factor  associated 
4 

The  parameter  correction  signal  Y^(t)  has  the  same  form  as  the  parameter 
drift  and  has  some  additional  phase  lag  compared  to  the  parameter  drift. 
This  correction  signal  will  be  fed  back  negatively  to  parameter  adjustment 
and  corrects  the  drifted  parameter  to  the  optimum  parameter  value. 

(b)  For  the  step  parameter  drift  case3  Y^(t)  equation  (3*23)  is  multiplied 
by  Yg(t)  equation  (3-24)  and  we  get 

Y_  (t )  =  A  DCB  C,sin(  IV  -  0  )  -  A,D.B  C-Cos^^  t  -  0O) 

7  0  0  1  1  p  I  o  1  I  z 
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(3-29) 


+  A  D  B  C,(X  -  X.)  sin2 (  U>,t  -  6  ) 

J  /  i  O  1  p 

In  Y^(t)  only  the  last  term  develops  a  D.  C.  component,  which  will  be 
integrated  and  drive  the  parameter  adjustment.  Other  terms  do  not  contri¬ 
bute  to  parameter  correction  because  they  are  comparatively  high  frequency 
sinusoids  and  and  are  very  small  due  to  the  bandpass  filter  attenua¬ 
tion  and  will  be  further  attenuated  by  the  smoothing  filter.  If  only  the 
last  term  of  Y^(t)  equation  (3-29)  is  expanded  it  becomes  as, 

-i—  A3D?B  CjCX-Xq)  -  4"A3D7B  C^X-Xq) 0052(^1  -  0  )  (3-30) 

Putting  back  the  original  value  for 

Y7d.c.(t>  -  iiVp‘)bc!<*  -  v  (>31) 

Thus  the  rate  at  which  the  parameter  drift  is  corrected,  is  proportional  to 
the  parameter  error  (X-X^),  the  perturbation  signal  amplitude  the  con¬ 
stant  a^  associated  with  the  IP  curve  (the  sersitivity  of  the  IP  variation 
to  changes  in  parameter  value)  and  adaptive  loop  gain.  This  open  loop 
analysis  is  not  exactly  correct  when  the  loop  is  closed,  but  it  shows  a 
possible  signal  development  through  the  adaptive  loop. 

(c)  Next,  the  effect  of  phase  difference  between  the  carrier  and  the  detec¬ 
tion  reference  signal  Y,(t),  in  the  demodulation,  is  considered.  To  see 

o 

more  clearly  the  process  of  demodulation  and  the  effect  of  the  phase  dif¬ 
ference  between  carrier  and  the  detection  reference  signal,  graphical  multi¬ 
plication  of  actual  signals  involved,  is  performed.  Taking  the  case  of 
step  parameter  drift,  the  actual  modulated  signal  in  Y^(t)  has  the  following 
forms  as  explained  in  3.2(b)  and  expressed  in  equations  (3-16)  and  (3-17), 


21 


A  |  X  -  X0  |  sinC&^t  -  6  )  when  X  -  XQ>  0 


(3-32) 


A  I  X  -  Xn  I  sin(  t  -  G  -  n)  when  X  -  0 

I  U  |  Ip  u 


(3-33) 


where  A  is  a  constant* 

The  carrier  in  Y^ft)  is  C^5(t)  =  sin(  CV^t  -  0  )  when  X~XQ  >  0  and  -C  ^(t)  = 
sin(<^>  t  -  0  -  it)  when  X-X.,  <*  0* 

The  detection  reference  signal  is  Y^(t)  =  B  sin(  CO  t  -  0  ) . 

If  G^  ^  0^  Y^(t)  leads  or  lags  the  carrier  C  ^(t)  by  some  angle.  The 

graphical  multiplication  is  performed  for  the  following  three  cases, 

(1)  Yg(t)  is  in  phase  with  the  carrier. 

(2)  Y, (t)  leads  the  carrier  by  30°. 

o 

(3)  Y,(t)  lags  the  carrier  by  30°. 

b 

The  coefficients,  of  each  term  is  a  constant*  therefore  only  a  sine  wave  of 

unit  maximum  amplitude  represents  each  wave, 

Fig.  5  shows  the  demodulated  output  by  dotted  line  for  the  above  three 

cases.  It  shows  clearly  the  phase  sensitive  nature  of  perturbation  signal; 

when  the  carrier  in  a  modulated  signal  is  C  j(t)  -  sin ( -  0^)*  the 

demodulated  output  is  all  positive  and  when  the  carrier  is  -C  -(t)  =  sin 

(  tO .  t  -  0  -  Jt)  the  demodulated  output  is  all  negative.  Thus  it  will  give 

1  P 

correct  direction  of  parameter  correction  in  accordance  with  the  actual 
parameter  error*  X-X^  >  0  or  X-X0  <  0. 

In  case  (1),  we  get  maximum  net  output  of  correction  signal  when  inte¬ 
grated  as  seen  in  Fig.  5(1).  For  the  case  (2)*  we  get  less  net  output  of 
correction  signal  when  integrated.  This  eventually  ends  up  with  less  adap¬ 
tive  loop  gain.  The  case  (3)  has  the  same  output  as  in  the  case  (2), 

If  we  consider  the  time  lag  of  starting  to  build  up  correction  signal 
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output  at  the  integrator,  the  case  (3)  has  the  largest  time  lag  of  the 

three  cases,  and  the  case  (1)  is  the  second  and  the  case  (2)  is  the  least* 

But  this  time  lag  is  negligible  because  the  carrier  frequency  is  very  high 

compared  to  the  modulating  frequency* 

From  this  analysis,  we  can  see  that  when  the  phase  difference  between 

Yc  and  Y.  is  exactly  90° ,  the  signal  output  will  be  the  lowest  and  when  the 
d  b 

phase  difference  is  more  than  90°  the  demodulated  signal  polarity  will  be 
reversed  and  the  magnitude  will  increase  as  the  phase  difference  approaches 
180° •  Therefore,  when  the  carrier  phase  lag  is  expected  to  be  closer  to  180°, 
it  is  better  to  remove  the  reference  signal  phase  lag  filter  and  change  the 
polarity  of  the  demodulated  output  to  have  correct  feedback 
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4.  The  Simulation  Study  of  a  Sinusoidal  Perturbation.  System 

4.1  IP.  variation  investigation 

To  investigate  the  actual  variation  of  the  IP.  with  respect  to  the 
parameter  drift  and  perturbation  of  the  plant s  the  following  model  and 
plant  were  simulated  by  an  analog  computer  with  the  adaptive  loop  opened. 
The  analog  computer  used  throught  the  investigation  is  the  Donne r  3100. 


M0del;  Gmod(s) 


K 

S(S  *  x0) 


Plant;  G  (s)  *  yyr - r 

p  S{S  +  x) 


The  pole  of  the  plant  is  a  variable  parameter  xs  and  the  pale  of  the  model 
is  fixed  at  3,  which  is  the  optimum  parameter  value  in  this  case  and  gives 
J  of  0.5.  The  pole  of  the  plant  can  vary  in  the  range  that  gives  J f  of 
plant  from  0.1  to  1.  with  K  ~  9.  The  block  diagram  of  the  simulation  is  as 
Fig.  6  and  the  analog  computer  circuit  is  shown  in  Fig.  7. 

(a)  The  input  to  the  plaint  and  model  Is  a  square  wave  of  0.1  cps. 

A  plant  parameter,  the  pole  of  the  plant  is  drifted  by  a  triangular  wave 
of  0.01  cps.  while  the  model  and  plant  are  operating.  The  output  of  model 
and  plant  are  as  Fig.  8(a)  and  (b)  respectively.  The  output  wave  shape 
shows  the  change  of  jf  of  the  plant  due  to  the  parameter  drift. 

(b)  The  error  output  is  recorded  as  Fig.  8{c)  and  the  pole  drift 
is  as  Fig.  8(d). 

Next,  the  input  to  plant  and  model  is  sine  v&ve  of  0.3  cps.  and  the 
error  output  is  also  recorded  as  Fig.  9{a).  The  Fig.  9{b)  is  the  pole  drift 
as  before  and  (c)  is  the  error  squared  output. 

The  above  results  shows  that  error  is  produced  only  when  the  output 
of  plant  is  in  transient  and  the  plant  response  is  net  the  same  as  the 
model  response  due  to  the  parameter  drift..  Also  it  shews  that  the  magnitude 
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Fig*  6  The  model  and  plant  set-up  for  the  investigation 
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Fig*  7  The  analog  computer  circuit  of  the  model  snd  plant  set-up  for  the 
investigation  of  IF. 


Pig.  8.  The  output  of  model,  plant  and  errorwith  plant  pole  drifted. 

(a)  model,  ■’ 

(b)  plant 

( c )  error 

(d)  plant  pole  drift 
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Pig.  9.  Error  and  error  squared  output  with  sine  wave  input  to  plant, 
fa)  error 

(b)  plant  pole  drift 
(_c)  error  squared 
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of  error  is  proportional  to  the  parameter  error,  the  magnitude  of  para¬ 
meter  drift. 

The  magnitude  of  error  is  bigger  when  the  parameter  drift  is  in  the 
direction  of  decreasing  plant  with  the  same  amount  of  parameter  drift. 
Therefore,  it  is  expected  that  the  magnitude  of  parameter  correction  signal 
will  be  greater  for  the  parameter  drift  to  a  decreasing  than  a  increasing 
J9  ,  with  the  same  amount  of:  drift, 

(c)  Error  output  is  observed  with  only  perturbing  the  pole  of 
plant  by  a  sinusoidal  perturbation  signal.  The  perturbation  frequency  is 
much  higher  than  a  drift  frequency;  here  a  sine  wave  of  5  cps  is  used. 

(1)  The  magnitude  of  error  output  due  to  perturbation  signal 
only.  The  perturbation  wave  appears  in  the  error  output,  but  the  magnitude 
is  much  reduced  due  to  the  severe  attenuation.  Five  cps  is  a  high  frequency 
for  this  plant.  Thus  the  parameter  perturbation  signal  does  not  appear  con¬ 
spicuously  in  the  plant  output  even  with  a  big  amplitude  perturbation  signal 
as  far  as  the  frequency  of  perturbation  is  high.  If  we  want  to  eliminate 
the  effect  of  perturbation  on  plant  output  completely,  the  parameter  of 
model  can  be  perturbed  instead  of  plant  parameter.  Fig,  10(a)  is  the  per¬ 
turbation  wave  appearing  on  the  output  of  the  plant.  The  perturbation 
waves  appear  conspicuously  on  the  output  of  the  plant  only  when  the  ratio  of 
the  amplitude  of  the  plant  input  to  the  amplitude  of  the  perturbation  signal 
is  small, 

(2)  Phase  shift  of  perturbation  signal  through  plant. 

The  perturbation  wave  appearing  on  the  output  of  plant  lags  the  original 
perturbation  signal  by  some  angle.  This  phase  lag  is  introduced  by  the 
plant  and  depends  on  the  characteristic  of  plant  and  signal  frequency. 

Fig.  10(c)  and  (d)  are  the  original  perturbation  signal  and  the  perturbation 
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(a)  Perturbation  wave 
appearing  on  plant 
output. 


(b)  Original  perturbation 
signal. 


(c)  Original  perturbation.  (e)  Input  to  plant 

(d)  perturbation  wave  on  output  (f )  Output  of  plant 

Pig.  10.  Perturbation  waves  on  the  output  of  plant. 


31 


wave  appearing  on  the  plant  output  respectively.  The  same  marks  on  the 
waves  are  the  corresponding  points  on  both  waves.  These  corresponding 
points  are  found  by  starting  and  stopping  the  signal  abruptly.  According 
to  the  measurement  of  phase  difference  between  the  two  waves,  perturbation 
wave  appearing  on  the  output  lags  the  original  signal  by  165°.  To  compare 
this  phase  lag  with  the  phase  lag  introduced  by  plant  to  a  signal  coming 
into  the  plant  through  the  plant  input,  the  same  signal  used  as  perturbation 
is  fed  into  the  plant.  Fig.  10(e)  and  (f)  are  the  input  to  plant  and  the 
output  of  plant.  The  same  marks  on  the  signals  are  the  corresponding  points 
on  both  waves  which  is  found  as  before.  According  to  the  measurement,  the 
output  waves  lags  input  by  165°  approximately,  which  is  the  same  phase  lag 
introduced  by  plant  to  the  parameter  perturbation  signal.  From  this  re¬ 
sult  it  is  concluded  that  the  phase  lag  introduced  by  the  plant  to  the 
parameter  perturbation  signal  fed  through  the  pole  of  plant  is  the  same 
phase  lag  introduced  by  plant  to  a  signal  of  same  frequency  fed  through  the 
plant  input. 
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4.2  The  Simplest  Adaptive  System  Simulation  (Adaptive  circuit  1) 

The  complete  adaptive  control  system  is  as  Fig.  1  and  the  block  dia¬ 
gram  is  as  Fig.  4.  To  investigate  the  function  of  each  component  and  the 
effect  of  change  in  components  on  the  over-all  system  operation,  the  simula¬ 
tion  is  started  with  the  simplest  combination  of  components  for  a  adaptive 
system.  Then  more  components  are  added  to  improve  the  system  operation 
until  the  complete  adaptive  system  is  formed.  We  designate  successive  cir¬ 
cuits  by  number  from  the  simplest  to  the  complete  circuit  in  simulation 
study.  The  first  circuit  is  designated  as  adaptive  circuit  1,  which  is  the 
simplest  and  is  shown  in  Fig.  11(a).  The  analog  computer  set-up  is  as  Fig, 
11  (b).  It  simulates  starting  from  the  error  output  point  in  the  complete 
system.  In  Section  3.2,  signal  analysis,  the  instant  deviation  of  parameter 
Y^(t)  is  expressed  as  equation  (3-4), 

Yj (t)  =  X(t)  +  Cp  -  XQ  (4-1) 

For  the  analysis  the  optimum  parameter  value  is  arbitrary  and  can  be 
set  as  zero.  Y^(t)  causes  an  error  output  and  the  error  output  has 

severe  attenuation  and  phase  shift  and  comes  out  as  equation  (3-6)  for  a 
sinusoidal  parameter  drift.  But  in  the  simulation  study,  ignore  the  at¬ 
tenuation  and  phase  shift  and  assume  the  error  output  Y2(t)  as, 

Y,(t)  =  X(t)  +  C  (4-2) 

^  P 

This  error  output  is  going  to  be  the  input  in  the  Adaptive  circuit  1  because 
the  circuit  simulates  from  error  output  point.  In  an  actual  system,  the 
parameter  correction  signal  corrects  plant  parameter  and  reduces  the  error 
output  to  zero  but  in  this  simulation  the  correction  signal  directly  reduces 
the  Y2(t),  which  is  equivalent  to  correcting  the  plant  parameter. 
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(a)  block  diagram 


Pig.  11.  Adaptive  circuit  1. 
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In  the  circuit,  and  are  a  squarer  and  a  detector,  both  of  which 

are  electronic  multipliers.  The  zero  D.  C.  level  of  electronic  multiplier 

is  not  constant,  therefore  a  pure  integrator  can  not  be  put  in  right  after 

the  multiplier  M«  because  the  integrator  integrates  the  drifting  D.  C.  and 

K4 

developes  erratic  correction  signal.  A  pseudo  integrator  - - is  used 

v  s  +  u .  1 

4 

for  the  simulation  of  the  pure  integrator  — - -  .  The  parameter  perturba- 

b 

tion  and  parameter  drift  signals  used  are  as  follows: 

0^  =  CjSin^jt  Cj  =  1  v,  ^  =  5  cps. 

X(t)  =  X.sin£c>  t  X,  =  8  v.  CU  =0.1  cps. 

I  m  1  m 

The  adaptive  loop  is  open  at  Yg  feedback  point  and  open  loop  observation  is 
made . 

(a)  The  wave  shapes  of  signals  at  every  point  of  loop  are  shown 
in  Fig.  12.  The  (a)  in  the  Fig.  is  the  error  signal  in  which  perturbation 
signal  is  superimposed  on  the  parameter  drift  signal.  (b)  is  the  output  of 
squarer  where  the  modulation  is  accomplished.  (c)  is  the  output  of  detector; 
the  demodulated  output  has  positive  or  negative  D.C.  component  in  accordance 
with  the  parameter  drift  direction.  This  is  integrated  and  comes  out  as  (d), 
which  is  the  parameter  correction  signal.  If  the  parameter  correction  sig¬ 
nal  is  compared  with  the  parameter  drift  signal  in  (a),  it  has  the  same  shape 
as  the  parameter  drift,  therefore,  when  the  correction  signal  is  negatively 
fed  back  it  will  correct  the  parameter  drift.  But  the  correction  signal  lags 
the  drift  signal  by  60°  approximately  as  seen  in  (a)  and  (d).  The  marked 
points  are  the  corresponding  point  on  both  waves.  This  phase  lag  is  caused 
by  the  integrator  as  seen  between  (c)  and  (d).  There  is  no  phase  lag  be¬ 
tween  Yj  and  Y^  as  seen  in  (a),  (b)  and  (c).  The  perturbation  signal  ap¬ 
pears  in  the  correction  signal  Yg,  which  will  be  eliminated  by  adding  a 
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(b)  Error  squared 
Y?(t) 


(a)  Error,  Y2(t) 


(c)  Demodulated 
Y?(t) 

-10 


(d)  Correction  signal  0 


Y9(t) 


CHART  NO.  RA  2921  32  BRUSH 


Pig.  12.  Wave  shapes  of  signals  in  Adaptive  circuit  1. 
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smoothing  filter, 

A  square  wave  also  can  be  used  as  a  parameter  perturbation  signal. 

The  same  results  as  obtained  by  sinusoidal  perturbation,  are  obtained 
using  square  wave  perturbation.  The  square  wave  perturbation  wave  ap¬ 
pearing  in  the  correction  signal  looks  like  a  sinusoidal  wave.  Fig.  12-1 
shows  the  results. 

(b)  When  the  frequency  of  perturbation  is  increased,  it  is  at¬ 

tenuated  more  and  it  does  not  appear  in  the  correction  signal;  in  this 
case  there  is  no  perturbation  signal  in  the  correction  signal  when  15  cps. 
is  used  for  The  amplitude  of  the  correction  signal  is  not  affected 

by  the  change  of  perturbation  signal  frequency.  This  means  the  perturba¬ 
tion  frequency  does  not  affect  the  loop  gain  in  this  simple  circuit,  how¬ 
ever,  in  an  actual  circuit  the  frequency  of  perturbation  will  affect  the 
amount  of  attenuation  due  to  circuit  components  and  affect  the  loop  gain. 

Fig.  13  shows  the  correction  signal  with  different  perturbation  frequencies. 

(c)  Next,  the  perturbation  frequency  is  fixed  at  5  cps.  and  the 
frequency  of  parameter  drift  is  varied.  When  the  drift  frequency  increases, 
the  amplitude  of  correction  signal  is  decreased.  Fig.  14  is  the  amplitude 
of  correction  signal  vs.  drift  frequency  graph.  The  decrease  of  amplitude 
of  correction  signal  is  due  to  the  attenuation  by  the  integrator.  This  re¬ 
sult  and  the  result  from  (b)  show  that  in  a  closed  loop  analysis  of  the 
adaptive  circuit,  we  need  to  deal  with  the  parameter  drift  signal  only. 

The  perturbation  signal  does  not  affect  the  magnitude  of  correction  signal 
even  though  it  comes  through  the  detector  and  rides  on  the  correction  signal. 
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(a)  Error, Yg(t)  0 
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Fig.  12-1  Wave  shapes  of  signals  in  adaptive  circuit  1  with  square 
wave  perturbation* 
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(a)  when  CO 


5 


CHART  NO.  RA  2921  32 


BRUSH  INST 


cps.  0 


10  sec* 


(c)  when  CO 
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Pig.  13.  Correction  signal  with  different  perturbation  frequency. 
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4*3  Adaptive  circuit  2 


From  the  simulation  study  of  adaptive  circuit  Ij,  it  is  found 
smoothing  filter  is  needed  to  eliminate  perturb&t io»n  signal  appearing  in 
correction  signal*  A  smoothing  filter  G^(s)  is  inserted  between  and 
the  integrator  in  adaptive  circuit  1*  The  circuit  becomes  as  Fig*  15(a) * 

The  analog  computer  circuit  is  as  Fig*  15(b)*  This  circuit  is  designated 
as  adaptive  circuit  2*  The  smoothing  filter  will  filter  out  completely  the' 
perturbation  signal  appearing  on  the  correct  ion  signal  but  also  attenuates 
the  correction  signal  and  gives  more  phase  lag  to  the  correction  signal* 

Too  much  phase  lag  of  the  correction  signal  tends  to  m&e  a  closed  loop 
system  unstable*  It  causes  time  lag  of  the  correcting  parameter  and  when 
the  parameter  change  is  abrupt  such  as  step  change^  excessive  overshoot  re¬ 
sults  on  the  correction  signal  and  sometimes  it  produces  a  sharp  undershoot 
before  producing  correction  signal  in  correct  direction*  The  undershoot  is 
a  strong  correction  signal  in  a  wrong  direction*  It  will  be  studied  in  more 
detail  later*  The  smoothing  filter  is  usually  chosen  as  a  compensator  for 
the  adaptive  loop  to  improve  step  response  cf  the  system* 


(a)  In  the  simulation£ 


filter  of  the  form. 


+  T2s 


is  used* 


The  *2* £  valtaes  used  are  0*2^  0*59  ls  2S  5S  10*  The  perturbation  signal 

and  parameter  drift  signal  used  in  4*23  are  used*  Taking  the  perturbation 
signal  present  in  correction  signal  as  noise^  correction  signal  to  noise  ratio 
is  measured  at  the  output  point  of  the  smoothing  filter  Fig*  15  is 

the  graph  of  correction  signal  to  noise  ratio  vs*  of  the  smoothing  filter* 
According  to  the  graph 9  'Z^  of  2  or  above  gives  more  than  30  :  1  signal  to 
noise  ratio*  The  perturbation  waves  remaining  in  tb*  correction  signal  even 
after  the  filter  ing9  will  fee  further  attenuated  by  the  int  egret  or  following* 
of  2  gives  corner  frequency  of  0.5  rad. /sec.  *rd  this  gives  enough 
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(a)  Block  diagram 


(b)  Computer  circuit 
Fig.  15.  Adaptive  Circuit  2. 
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Signal 

to 

noise 

ratio 


Fig.  16  The  signal  to  noise  ratio  vs.  filter  time  constant 
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at tenuat ion  to  perturbat ion  waves  of  5  cps.  present  in  correction  signal  of 
this  system.  Thus  it  is  concluded  that  an  acceptable  ratio  of  corner  fre¬ 
quency  of  smoothing  filter  to  perturbation  frequency  is  1  :  60  or  greater. 

But  the  attenuation  of  perturbation  waves  in  the  correction  signal  is  not  so 
important  as  the  stability  of  the  system,  therefore,  this  filter  is  designed 
as  an  appropriate  compensator  for  the  adaptive  loop  and  it  naturally  acts 
as  a  smoothing  filter. 

(b)  Effect  of  adaptive  loop  gain  on  parameter  drift  correction 

(1)  Sinusoidal  parameter  drift 

The  parameter  is  drifted  sinusoidally.  When  the  adap¬ 
tive  loop  is  open  the  parameter  drifts  sinusoidally  and  as  soon  as  the 
adaptive  loop  is  connected  the  drifting  parameter  is  corrected  to  the  optimum 
value.  The  loop  gain  affects  the  amount  of  steady  state  error  as  expected. 

The  necessary  loop  gain  is  decided  by  the  required  steady  state  accuracy. 

Fig.  17(a)  shows  a  parameter  drift  when  the  adaptive  loop  is  open  and  then 
corrected  parameter  condition  when  the  loop  is  closed.  (b)  is  the  corres¬ 
ponding  parameter  correction  signal.  The  loop  gain  is  a  little  low  and  the 
steady  state  error  is  comparatively  large.  To  reduce  the  steady  state  error, 
loop  gain  is  increased  and  the  parameter  drift  is  almost  completely  correct¬ 
ed  as  seen  in  Fig.  17(c).  With  this  high  loop  gain,  the  loop  was  closed 
suddenly  when  the  parameter  drift  was  a  maximum,  the  parameter  drift  was 
corrected  to  the  optimum  value  with  an  excessive  oscillation  before  settling 
down  to  steady  state;  Fig.  18(a)  shows  the  parameter  drift  being  corrected 
and  (b)  is  the  corresponding  correction  signal  wave  shape. 

(2)  Triangular  Parameter  Drift 

Fig.  18(c)  and  (d)  is  the  triangular  parameter  drift  being 
corrected  and  its  correction  signal.  The  loop  gain  is  a  little  low  and  the 
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10 

(c)  0 
-10 


Pig.  17.  Parameter  drift  correction 

fa)  Parameter  drift  being  corrected  with  low  loop  gain 

(b)  Parameter  correction  signal 

(c)  Parameter  drift  being  corrected  with  a! high' loop- gain 
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(c) 


(4) 


Pig.  18.  Parameter  drift  correction 

(a)  Sinusoidal  parameter  drift  being  corrected  with  high  gain 

(b)  Correction  signal 

(c)  Triangular  parameter  drift  being  corrected 

(d)  Correction  signal 
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correction  is  not  perfect. 

(3)  Square  wave  parameter  drift 

The  parameter  is  drifted  by  a  square  wave.  Fig.  19(a) 
shows  the  square  wave  parameter  drift  being  corrected  and  (b)  is  the  corres¬ 
ponding  correcting  signal,  in  which  steady  state  error  is  small.  The  system 
is  stable.  When  the  loop  gain  is  too  high,  the  system  response  becomes 
oscillatory  and  tends  to  be  unstable  but  the  system  will  never  be  unstable 
because  overall  loop  is  second  order  system  in  this  set-up.  Fig.  19(c) 
is  the  parameter  being  corrected  and  (d)  is  the  corresponding  correction 
signal  when  the  loop  gain  is  very  high. 

(c)  The  effect  of  perturbation  frequency  on  adaptive  loop  step 

response . 


The  most  extreme  case  of  parameter  drift  is  step  change  of 
parameter  due  to  an  abrupt  environmental  change.  Therefore,  only  the  effect 
of  perturbation  frequency  on  adaptive  loop  step  response  is  investigated. 

t 

The  plant  parameter  is  drifted  by  square  wave  of  0.01  cps.  For  dif¬ 
ferent  frequencies  of  perturbation,  maximum  oveishoot  and  rise  time  are 
measured.  The  perturbation  frequencies  are  varied  from  1  cps.  to  100  cps. 

(1)  Steady  state  magnitude  of  correction  signal 

In  this  adaptive  circuit  there  is  no  attenuating  com¬ 
ponent  for  the  perturbation  signal  between  perturbation  signal  input  and  the 
demodulation,  therefore,  the  perturbation  frequency  does  not  affect  the  loop 
gain  and  the  steady  state  magnitude  of  correction  signal  is  constant  with 
different  perturbation  frequency  as  seen  in  Fig.  13. 

(2)  Maximum  overshoot  and  undershoot 

For  a  very  low  frequency  of  perturbation,  the  correction 
signal  has  extremely  irregular  wave  shapes  of  square  waves  as  shown  in  Fig. 
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Pig.  19.  Parameter  drift  correction 

(a)  Square  wave  parameter  drift  being  corrected 

(b)  Corresponding  correction  signal 

(c)  Square  wa-se  parameter  drift  being  corrected  with  very  high 
loop  gain 

(d)  Corresponding  correction  signal 
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(b)  corresponding  correction  signal. 


20(b),  The  perturbation  frequency  is  2  cps.  The  correction  signal  has 
various  amounts  of  overshoots  and  undershoots.  Here,  undershoot  is  de¬ 
fined  as  a  sharp  and  strong  correction  signal  to  a  wrong  direction.  The 
magnitude  of  overshoot  and  undershoot  depends  on  the  point  of  perturba¬ 
tion  signal  wave  at  which  an  abrupt  change  of  parameter  occurs.  The  ir¬ 
regularities  of  correction  signal  wave  shapes  becomes  greater  at  lower 
perturbation  frequency  and  with  perturbation  frequency  below  1  cps.  the 
correction  signal  is  too  irregular  to  have  a  sensible  correction.  The 
magnitude  of  overshoot  and  undershoot  also  decreases  as  the  perturbation 
frequency  increases  and  above  5  cps,  of  perturbation  frequency  there  is 
no  overshoot  nor  undershoot.  Table  4,1  lists  the  measured  ratio  of  maxi¬ 
mum  overshoot  and  undershoot  to  steady  state  correction  signal  amplitude  at 
lower  perturbation  frequencies, 

C0^  (cps ,  )  Max,  overshoot  Max.  undershoot 


1 

1.63 

4.25 

1.5 

1.28 

2.06 

2 

1.19 

1.81 

3 

no  overshoot 

1.38 

4 

tt 

1.18 

5 

II 

1.03 

Table  4,1  Ratio  of  Maximum  overshoot  and  undershoot  to  the  amplitude 
of  steady  state  correction  signal. 

It  is  noted  from  Table  4,1  that  the  maximum  undershoot  ratio  is  much  great¬ 
er  than  overshoot  ratio.  The  sharp  undershoot  should  be  avoided,  which 
makes  the  system  unstable  and  leads  to  erratic  operation.  Therefore,  the 
perturbation  frequency  should  be  higher  than  5  cps,  in  this  adaptive  control 


system. 
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(3)  Rise  time 

Rise  time  vs,  perturbation  frequency  is  plotted  as  Fig. 
21.  The  rise  time  increases  sharply  from  CO^  =  1  cps,  to  aJ  «  4  cps. 

Above  5  cps,  rise  time  increases  very  little  and  stays  almost  constant. 

Fast  rise  is  desirable  and  it  is  obtained  at  lower  Co  but  it  accompanies 
undesirable  undershoot  as  seen  in  Table  4.1.  In  the  following  simulation 
studies,  CO^  -  5  cps,  is  used  as  perturbation  frequency  throughout,  which  is 
the  lowest  perturbation  frequency  that  does  not  accompany  any  overshoot  or 
undershoot  in  correction  signal. 
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4.4  Adaptive  Circuit  3 

A  band  pass  filter  G^(s)  was  added  between  the  squarer  and  the  de¬ 
tector  in  adaptive  circuit  2,  Fig.  15.  The  circuit  becomes  as  shown  in 
Fig.  22.  The  band  pass  filter  attenuates  all  unnecessary  frequency  compon¬ 
ents  considerably  in  the  signal  coming  from  squarer  and  passes  the  components 
which  contain  information  without  attenuation.  The  band  pass  filter  is  es¬ 
sential  in  a  multi-loop  adaptive  control  system  where  different  perturbation 
frequencies  carry  information  for  different  parameters  being  controlled  through 
different  loops.  This  filter  also  affects  the  stability  of  the  adaptive  loop. 
The  requirement  of  the  band  pass  filter  is  explained  in  Section  3.3.  To  meet 
the  requirement  a  filter  of  following  form  is  used, 


a?  s 


^i(s)  =  o  2 

S  +  2f(t>S+£0{ 

where  =  perturbation  frequency  (center  frequency) 
^  =  damping  ratio 

The  half  power  frequency  of  this  filter  is  calculated  as, 


(4-3) 


(4-4) 


The  band  width  46>b.w.  =  2 J  CQ ^ 

To  decide  the  value  which  gives  best  step  response  to  the  system,  the  T* 

is  varied  and  the  step  response  of  the  system  is  observed.  jJ*  of  0.5  or 
0.6  gives  best  response. 

(a)  First,  only  the  band  pass  filter  is  simulated.  ^  of  0.5  is  taken 
and  the  filter  becomes  as, 

4y2s 

l 


G1(s)  = 


s2  +  S  +  &>2 


(4-5) 
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G4(s)  G2(s) 

Pig.  22.  Adaptive  circuit  3 


Fig.  23.  Analog  computer  simulation  of  band  pass  filter  G^s) 
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The  analog  computer  simulation  of  the  band  pass  filter  is  shown  in  Fig,  23, 
The  frequency  response  and  phase  angle  characteristic  is  plotted  as  Fig,  24. 
The'  calculated  half-power  frequency  and  band  width  and  the  actual  half¬ 
power  frequency  and  band  width  are  as  below, 


(rad/sec) 

&7h2(  rad/sec) 

CO ,  „  (rad/sec) 

calculated 

15.7 

47.1 

31.4 

actual 

19.2 

49.5 

30.4 

Using  the  result  of  signal  analysis  in  Section  3.3(a),  following  signal 


components  are  expected  to  be  present  in  the  signal  from  the  squarer  of 
adaptive  circuit  3, 

D.C 

2  CO  =  1.256  rad/sec 
m 

2  U)  c  62.8 

CO.  -  U)  s  30.144  rad/sec 
1  m 

UJn  +  CO  =  32.656  " 

1  m 

Therefore,  this  filter  will  suppresses  D.C.,  2  CO^,  and  2  Ci)^  and  passes 
(/).- (jQ  ,  and  (jO  +  44?  components  without  attenuation. 


(b)  The  adaptive  circuit  3  simulation  was  obtained  by  connecting 

the  band  pass  filter  simulation,  Fig.  24  between  and  M2  in  adaptive 

circuit  2,  Fig.  15.  The  performance  of  the  adaptive  loop  is  much  improved 

by  the  addition  of  the  band  pass  filter.  With  proper  loop  gain,  a  square 

wave  parameter  drift  is  corrected  perfectly  as  shown  in  Fig.  25,  where  (a) 

is  the  square  wave  parameter  drift  being  corrected  and  (b)  is  its  correction 

signal.  The  H  is  1.08  and  the  rise  time  is  0.7  sec.  Even  with  excessive 
pt 

loop  gain,  the  performance  is  less  oscillatory  compared  to  the  same  case  in 
adaptive  circuit  2,  in  4.3(b)(3).  Fig.  25(c)  is  the  parameter  drift  being 
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tTI±l 


Pig.  25.  Parameter  drift  correction 

(a)  Square  wave  parameter  drift  being  corrected 

(b)  Corresponding  correction  signal 

(c ^Square  wave  parameter  drift  being  corrected  with  excessive 
loop  gain 

(d)  Corresponding  correction  signal 
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corrected  with  excessive  loop  gain  and  (d)  is  its  correction  signal* 

Compare  these  figures  with  Fig,  19(c)  and  (d)  which  are  the  same  case  in 
adaptive  circuit  2,  Fig*  25(c)  shows  that  the  parameter  drift  is  over¬ 
corrected  because  of  excessive  loop  gain  but  the  response  is  less  oscil¬ 
latory,  The  comparison  shows  a  definite  improvement  in  adaptive  loop  per¬ 
formance  by  adding  band  pass  filter.  This  filter  is  acting  as  a  compensator 
in  this  circuit, 

(c)  A  comparison  of  demodulated  outputs  is  made  between  following  two 
cases  in  sinusoidal  parameter  drift. 

(1)  Modulated  wave  is  not  filtered  by  band  pass  filter,  (adaptive 
circuit  2). 

(2)  Modulated  wave  is  filtered  by  band  pass  filter,  (adaptive 
circuit  3)* 

Fig*  26(a)  is  for  the  case  1.  Refer  to  the  Fig*  22  adaptive  circuit  3 
for  the  designation  of  signals  at  each  point*  For  the  case  1  the  modulated 
signal  goes  into  the  detector  without  filtering  by  band  pass  filter. 

The  demodulated  output  is  Y^.  The  modulated  signal  contains  many  un¬ 
necessary  signal  components  which  do  not  contain  any  intelligence  as  ex¬ 
plained  in  section  3,2,  Accordingly,  the  demodulated  output  contains  un¬ 
necessary  signal  components.  On  the  other  hand,  in  the  case  2,  Fig.  26(b), 
the  Y^  passes  the  band  pass  filter  which  filters  out  unnecessary  signal 
components  as  explained  in  section  3.3.  Thus,  the  signal  Y^  contains  only 
signals  which  contain  intelligence.  Note  the  wave  shapes  after  filtering. 
This  signal  is  demodulated  and  comes  out  as  Y^  which  has  clear. cut  polarity 
in  accordance  with  the  parameter  drift  direction.  This  indicates  the  Y7 
for  the  case  2  does  not  contain  any  unnecessary  signal  as  in  the  case  1. 

This  verifies  the  signal  analysis  in  section  3,3  and  3,4. 
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(a)  (b) 

Pig.  26.  Comparison  of  demodulated  output  with  modulated  signal  filtered 
and  unfiltered, 

(a)  Unfiltered  (b)  Filtered 
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4,5  Adaptive  circuit  4 

The  measurement  lag  filter  G  (s)  is  connected  between  the  parameter 

m 

perturbation  and  drift  input  and  squarer.  The  detection  reference  phase 
lag  filter  G^(s)  is  connected  between  reference  signal  input  and  the  de¬ 
tector.  The  Gm(s)  represents  the  phase  lag  in  measuring  error.  G^(s) 
chosen  to  give  the  same  phase  lag  to  the  detection  reference  signal  at 

as  G  (s)  gives  to  the  perturbation  signal  (the  carrier)  at  CV  ,  thus  the 
m  1 

carrier  and  the  detection  reference  signal  are  in  phase  with  each  other  at 
the  detector  The  band  pass  filter  G  (s)  does  not  introduce  an>  phase  lag 
to  the  carrier  at  60 and  the  total  phase  lag  introduced  to  carrier  be¬ 
tween  Y.  and  YcS  is  the  only  phase  lag  due  to  G  (s).  The  circuit  becomes 
1  D  m  i 

as  Fig.  27  and  analog  computer  simulation  is  as  Fig.  28.  This  circuit  is 

designated  as  adaptive  circuit  4.  The  filters  G  (s)  and  G«( s)  take  follow- 

m  3 

ing  forms, 


Gm(s)  "  1  +  r  S  ’ 


G3(S)  1  +  f.S 


(a)  The  process  of  demodulation  is  explained  and  analyzed  in  section 

3,4.  It  is  verified  in  this  simulation  study.  First,  ^  ~  ~  7 — 

is  used,  which  gives  45°  phase  lag  at  (O  =  31.4  rad/sec/ (5cps .) .  Thus 

the  carrier  signal  in  Yc  and  the  detection  reference  signal  Y.  are  in  phase 

j  o 

at  the  detector.  The  modulated  signal  Yc  is  multiplied  by  the  Y.  signal 
and  the  demodulated  output  is  produced.  The  demodulated  output  has  one 
polarity  or  the  other  in  accordance  with  the  parameter  drift  direction  with 
respect  to  the  optimum  parameter  value.  Fig.  29(a)  is  the  modulated  signal 
Y^  and  (b)  is  the  demodulated  output  Y7.  The  demodulated  output  has  exactly 

j  / 

the  same  wave  shape  as  obtained  by  graphical  multiplication  of  Yc  and  Y-  in 

o  o 


60 


O 
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Fig*  27  Adaptive  circuit 
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Fig*  28  Analog  computer  circuit  for  the  adaptive  circuit 


(b)  demodulated  output, 


section  3.4  Fig.  5.  This  signal  is  integrated  and  develops  a  correction 
signal  of  one  polarity  or  the  other.  When  the  polarity  of  detection  re¬ 
ference  signal  is  changed,  the  demodulated  output  also  changes  its 
polarity.  These  results  confirm  the  explanation  of  the  demodulation  in 
section  3.4. 

The  modulated  signal  and  demodulated  output  wave  shapes  for  square 
wave  parameter  drift  are  shown  in  Fig.  30.  The  (a)  in  this  figure  is  para¬ 
meter  drift,  (b)  is  the  output  of  squarer, (c)  is  the  modulated  signal  com¬ 
ing  out  of  the  band  pass  filter  and  (d)  is  the  demodulated  output.  The 
demodulated  output  has  one  polarity  or  the  other  depending  on  the  direction 
of  parameter  drift. 

(b)  Next,  the  effect  of  phase  difference  between  Y  and  Y,  on  the 

0  6 

demodulated  output  is  observed.  This  effect  is  explained  in  section  3.4 
and  the  result  of  graphical  multiplication  of  two  signals  which  are  out  of 
phase  with  each  other  by  some  angle  is  shown  in  Fig.  5.  The  response  of 
the  system  is  observed  for  the  following  three  cases. 

1.  Y,  and  Yc  are  in  phase. 

2.  Y.  leads  Yc  by  60°. 

o  D 

3.  Y^  lags  Y^  by  60°. 

For  the  case  one,  the  parameter  drift  being  corrected  and  parameter  cor¬ 
rection  signal  wave  shapes  are  as  Fig.  31(a)  and  (b)  respectively;  the  M 
is  1.046  and  rise  time  is  0.9  sec.  For  the  case  2  and  3,  the  results  are 
the  same  and  the  result  of  the  case  2  is  shown  in  Fig.  31(c)  and  (d).  The 
response  is  overdamped  and  the  rise  time  is  about  2  sec.  This  shows  that 
the  overall  loop  gain  is  reduced  due  to  the  phase  difference  between  the 
Yc  and  Y, .  This  effect  is  predicted  in  the  signal  analysis  in  section  3.4. 

The  phase  difference  between  Yc  and  Y,  does  not  affect  the  system  seriously 

o  o 

and  merely  reduces  the  loop  gain  which  can  be  easily  restored. 
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IVIEAITS  DIVISION  OF  CUVITE  CORPORATION  CLEV 


(b)  Output  of 
squarer 


(c) 


(d) 


time 


2 

0 

Demodulated 

output 

“2 


Modulated 
signal  after 
filtering 


Pig.  30.  Modulation  and  demodulation  of  square  wave  parameter  drift. 
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(v) 

5 

(a)  Parameter  0 
drift  being 
corrected  -5 


(c) 


(d) 


Pig.  31.  Effect  of  phase  difference  between  Yg  and  Yg 
(a),(b)  Yg  and  Yg  are  in  phase. 

(c),(d)  Yg  leads  or  lags  Yg  by  60° . 
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4.6  Complete  adaptive  control  circuit  (adaptive  circuit  5) 

Finally,  a  complete  adaptive  control  circuit  was  simulated  on  an 

analog  computer.  The  block  diagram  of  the  complete  system  is  as  Fig.  32. 

This  circuit  will  be  designated  as  adaptive  circuit  5,  The  analog  computer 

simulation  is  achieved  by  merely  replacing  the  G  (s)  in  the  adaptive  circuit 

m 

4  simulation.  Fig.  28  with  the  model  and  plant  simulation  circuit  in  Fig.  7 
and  the  correction  signal  is  negatively  fed  back  to  the  summer,  output  of 
which  goes  to  the  multiplier  M^.  The  reference  signal  phase  lag  filter 
G^(s)  is  not  used  because  the  carrier  is  expected  to  have  phase  shift  of 
160°  at  the  detector,  which  is  nearer  to  180°.  The  input  to  the  plant  and 
model,  is  R(t)  =  7sinl.884t  (f  -  0.3  cps.). 

(a)  The  smoothing  filter  ^2(3)  is  acting  as  a  compensator.  To  see  the 
effect  of  the  ^2(3)  on  the  system  response,  first  it  was  removed  and  the 
response  of  the  system  was  observed.  The  adaptive  loop  was  opened  at  the 
output  of  the  integrator  G^(s).  The  parameter  correction  signal  developed 
with  the  loop  open,  is  recorded  as  Fig.  33(b).  Fig.  33(a)  is  the  square 
wave  parameter  drift  signal.  To  close  the  adaptive  loop,  the  polarity  of 
the  correction  signal  should  be  reversed  to  have  negative  feedback. 

When  the  loop  is  closed,  this  system  becomes  unstable,  so  the  loop  gain 
was  reduced  to  a  very  small  value  such  that  the  system  barely  maintains  a 
stable  condition.  The  correction  signal  becomes  as  shown  in  Fig.  33(c). 

The  magnitude  of  the  correction  signal  becomes  very  small  and  the  signal 
shape  is  drastically  changed  compared  to  the  correction  signal  shape  with 
the  loop  open.  That  is  why  the  open  loop  signal  analysis  is  not  exactly 
valid  when  the  loop  is  closed  as  mentioned  in  section  3.  Fig.  33(d)  is 
the  plant  output  and  (e)  is  the  model  output.  Because  of  too  small  loop 
gain  and  accordingly  too  small  correction  signal,  a  slight  correction  of 
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32  The  complete  adaptive  control  circuit  (Adaptive  circuit  5) 


ED  IN  U  S  A 


Pig.  33  Signal  shapes  in  adaptive  circuit  5. 

(a)  parameter  drift,  (b)  Correction  signal  (open  loop) 
(c)  correction  signal  (closed  loop).\  (d)  plant  output 
(e)  model  output. 

l 
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parameter  drift  is  achieved.  Compare  the  figures  (d)  and  (e).  When  the 
plant  parameter  changes  abruptly  the  plant  output  has  sharp  overshoot  due  to 
instability  of  the  system, 

(b)  A  lag  filter  of  the  form  G9(s)  ^ -  is  inserted  be- 

Z  1  4*  l  2s 

tween  the  detector  and  the  integrator.  Three  different  values  of  2, 

5,  and  10  are  used  and  the  loop  gain  is  adjusted  to  give  stable  operation. 
With  different  values  of  'Z'^y  the  correction  signal  wave  shapes  are  differ¬ 
ent  but  in  all  cases  the  critical  loop  gain  is  too  low.  Thus  the  para¬ 
meter  drift  can  not  be  corrected  completely  in  this  circuit.  Fig*  34  is 
the  parameter  drift  being  corrected  and  corresponding  correction  signal 


for  the  three  cases;  Fig  34(a)  is  for  -  2,  (b)  is  for  ^  -  10,  (c)  is 
for  ^2  “  5,  An  effective  compensator  must  be  found  to  achieve  a  re¬ 
quired  performance  of  the  system. 
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(a) 


Fig*  34  Parameter  drift  and  correction  signal 

(a)?2  =  2  (b)  r2  =  10  (c)  r2  =  5 


\ 
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5.  Design  of  a  Sinusoidal  Perturbation  Adaptive  Control  System 

A  design  of  a  sinusoidal  perturbation  system  is  attempted  based  on  the 
simulation  study  and  using  a  linearization  technique.  The  design  procedure 
and  considerations  are  virtually  the  same  for  this  system  as  general  control 
systems  except  some  special  features  such  as  modulation  and  demodulation  of 
error  signal  using  perturbation  signal.  Finding  the  approximate  equivalent 
transfer  function  of  circuit  components  and  applying  linearization  technique, 
a  preliminary  design  is  made  on  analytical  basis  and  the  preliminary  design 
is  simulated  on  analog  computer  and  the  effectiveness  of  the  analytical 
method  is  checked.  The  final  refinement  of  design  should  be  done  on  experi¬ 
mental  basis, 

5.1  Specification 

There  are  many  ways  to  specify  the  performance  criteria.  Those  speci¬ 
fications  differ  depending  on  the  application 3  but  at  least  the  steady  state 
requirement  and  dynamic  or  transient  requirement  should  be  specified.  One 
special  specification  for  this  system  is  hunting  loss5  which  is  defined  as 
the  maximum  allowable  deviation  of  parameter  caused  by  sinusoidal  perturba¬ 
tion  at  the  optimum  parameter  value.  The  choice  of  the  hunting  loss  depends 
on  the  requirement  of  the  system  and  the  characteristic  of  the  IP  curve;  a 
small  hunting  loss  gives  less  distrubance  on  the  output  of  the  plant  but  may 
make  the  system  lock  into  small  local  minimum  in  IP  curve  and  sometimes 
does  not  give  any  information  about  the  actual  parameter  drift  because  of 
attenuation. 

5.2  Selection  of  Perturbation  frequency 

The  perturbation  frequency  CO ,  should  be  much  higher  than  the  expected 
maximum  frequency  of  the  parameter  drift  UJ^  this  is  necessary  in  any 
modulation  system.  The  investigation  of  the  effect  of  perturbation  frequency 
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on  the  system  response  in  4.3(b)  shows  that  a  perturbation  frequency  of 
higher  than  5  cps,  does  not  make  much  difference  on  the  system  response. 

At  an  extremely  low  frequency  of  perturbation,  the  system  operation  tends 
to  become  erratic  due  to  sharp  overshoot  and  undershoot  even  though  the 
rise  time  is  very  short.  If  the  perturbation  frequency  is  too  high,  the 
attenuation  through  the  system  components  is  severe  and  the  higher  gain  of 
perturbation  is  required. 

5.3  Linearization  of  the  System 

To  investigate  the  stability  and  response  of  the  system  analytically, 
an  equivalent  signal  transfer  function  of  the  adaptive  loop  is  established,, 
The  representative  circuit  block  diagram  is  derived  in  section  3.1  and 
shown  in  Fig.  4.  The  following  linearization  technique  is  applied  to 
linearize  the  adaptive  loop.  In  the  adaptive  loop,  a  D.C.  parameter  correc¬ 
tion  signal  is  proportional  to  the  plant  parameter  error  X  -  X^,  as  shown  in 
equation  (3-31)  section  3.4(b).  This  indicates  a  possibility  of  derivation 
of  an  equivalent  linear  loop.  In  applying  the  linearizing  technique  to  any 
specific  system,  only  the  variation  of  the  IP  curve  about  the  optimum  is 
important.  Therefore  in  the  general  IP  curve  equation,  set  X^  =  0  and  a^  = 

0.  Referring  to  the  Fig.  27  adaptive  circuit  4,  which  is  a  complete  block 

•I 

diagram  of  the  system,  the  portion  between  error  output  and  the  detector  is 
a  carrier  system.  To  linearize  the  squarer  and  the  detector,  the  describing 
function  or  a  first  term  of  Taylor  series  approximation  is  used  and  to  get 
an  equivalent  signal  transfer  function  for  band  pass  filter  G^(s)  frequency 
transformation  technique  is  used,  (Ref.  1) 

The  perturbation  signal  does  not  get  into  the  picture  when  a  closed  loop 
analysis  is  made  as  verified  in  simulation  study  4„20  Therefore,  only  the 
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modulating  signal  (parameter  drift  signal)  is  acting  in  the  closed  loop, 

(a)  The  parameter  drift  signal  acts  directly  on  the  G  (s)  and  no 

m 


linearization  is  necessary  for  G  (s). 

IQ 

(b)  Linearization  of 

Standard  describing  function  technique  is  used.  The  input  to  is  the 
error  output  Y ^ (t)  and  can  be  rewritten  from  equation  (3-7)  as, 


Y0(t)  =  A  C,sin(  frit  -  0  )  +  A  X-Sin(  CO  t  -  0  ) 
pi  x  1  p7  ml  m  nr 


(5-1) 


The  output  from  M2  is  Y^(t)  equation  (3-10),  In  this  equation  only  the 
last  term  is  the  carrier  modulated  by  the  intelligence  signal  and  all 
other  terms  does  not  contain  any  intelligence  and  will  be  suppressed  by 
the  band  pass  filter.  Therefore,  as  far  as  the  modulating  signal  is 
concerned,  only  the  last  term  of  the  Y^(t)  is  the  output  from  Mj,  The 
last  term  of  the  Y^(t)  is  written  from  equation  (3-10), 

2a.  A  AC  -  X1  s  in  (  CO  t  -  9  )sin  (  to  t  -  9 )  (5-2) 

impll  I  p  m  m 

The  equivalent  multiplier  gain  is  the  ratio  of  the  peak  amplitudes  of 

signal  input  and  output  terms.  Thus,  the  peak  amplitude  of  (5-2)  is 

divided  by  the  peak  amplitudeof  the  modulating  signal  term  in  (5-1)  and 

one  obtains 


G  -  =  2a  A  C. 
Ml  1  P  1 


(5-3) 


This  same  gain  also  can  be  obtained  by  a  Taylor  series  approximation. 

(c)  Linearized  equivalent  transfer  function  of  band  pass  filter  G^(s). 


Gj(s)  = 


2 

U)  S 


S2  +  +  CO* 
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The  steady  state  phase  shift  due  to  G^(s)  at  S  =  j  is  0. 

The  linearization  is  based  on  the  following  assumptions; 

1.  Phase  shift  curve  is  linear  at  S  =  J  OO  . 

2.  Attenuation  at  S  =  J(A>  +  CO  )  is  almost  the  same  as  at  S  =  1  CO,. 

j  >  j  —  J  1 

Considering  only  the  modulating  signal  component,  the  input  to  the  band 
pass  filter  is  the  signsl  (5-2)  and  the  output  is  the  last  term  signal  in 
Y^(t)  equation  (3-21).  The  output  is  written  as, 


2a  A  A  C.X  A  sin(  00  y t  -  9  )sin(  U)  t  -  0  -AG)  (5-4) 

1  m  p  1  1  3  1  p  m  m 

The  input  (5-2)  is  compared  with  the  output  (5-4),  the  modulating  signal 

is  attenuated  by  the  factor  of  A^  and  phase  lag  of  A  0  is  introduced  by 

the  G^(s).  An  equivalent  transfer  function  may  now  be  defined  such  as  to 

approximate  the  phase  characteristics  of  the  original  network  or, 


^1  equiv.^S^  1  +  ‘fr. 


(5-5) 


1  eq.’ 


In  this  case,  A  is  A^  which  is  the  attenuation  factor  of  the  band  pass 

filter  at  CO '  and  'Z,  is  determined  by  choosing  a  representative 

1  1  eq. 

value  oiOO  and  is  calculated  such  that  G,  .  (s)  gives  the  same 

m  1  eq.  lequiv. 

phase  lag  at  00  as  a  G,  which  results  in  the  original  system. 

(d)  Equivalent  transfer  function  for  M^. 

The  detection  reference  signal  input  to  Mj  is 


Y6(t)  =  BC1sin(^1t  -  0p)  (5-6) 

and  the  modulated  signal  input  to  M2  is  the  signal  term  (5-4). 

The  output  of  M,  is  the  Y  .  (t)  equation  (3-27)  and  is  rewritten  as, 

Y7d.c.(t)  =  *2  (2aiAmApClXlA3)BClSin(^mt  “  0m  *  A6)  (5-7) 

The  ratio  of  the  peak  amplitude  of  input  (5-4)  to  output  (5-7)  for  M2  is 
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(5-8) 


The  components  after  the  ^  in  the  adaptive  loop  act  upon  the  signal  terms 
directly,  so  no  further  transforms  are  required* 

5.4  Stability  investigation  of  uncompensated  system  and  compensation 

After  linearizing  the  required  components,  the  block  diagram  of  the 
system  is  formed  with  only  the  fixed  components*  All  components  are  fixed 
except  the  smoothing  filter  G2(s)  which  is  acting  also  as  a  compensator. 

The  block  diagram  is  reduced  to  a  convenient  form  and  the  overall  loop 
transfer  function  is  obtained.  The  gain  requirement  is  decided  on  the  basi 
of  static  and  dynamic  accuracy  and  the  gain  is  set  at  the  lowest  acceptable 
value.  Using  the  preferred  standard  method  of  analysis,  the  uncompensated 
system  is  analyzed.  The  minimum  gain  permissible  is  that  which  provides 
acceptable  steady  state  error.  In  most  uncompensated  systems  this  gain 
exceeds  the  critical  gain  of  the  system  and  the  system  becomes  unstable. 
Proper  compensation  is  necessary  to  make  the  system  stable  maintaining  the 
required  gain.  Many  different  methods  of  compensation  are  available  depend 
ing  on  situations,  but  whenever  possible,  a  phase  lag  compensation  is  recom 
mended  because  the  compensator  is  also  acting  as  a  smoothing  filter.  The 
required  band  width  is  narrow  in  this  system  because  the  parameter  drift 
frequency  is  generally  very  low  and  it  is  effective  in  suppressing  noise. 
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5,5  Example 

Some  adaptive  control  systems  are  designed  analytically  and  simulat¬ 
ed  on  computer.  The  actual  performance  of  the  system  is  compared  with  the 
one  predicted  by  analytical  method.  Because  of  difficulty  in  simulation 
on  an  analog  computer,  the  first  example  uses  mild  compensation  and  was 
simulated  on  analog  computer.  In  the  last  example  an  extreme  case  of 
design  is  made  and  it  is  simulated  on  digital  computer, 

(a)  Example  1 

Design  an  adaptive  control  loop  to  maintain  optimum  for  the 
following  plant  in  which  the  pole  is  drifting  with  environment  change. 


Cp(8)  . 


81 

S(S  +  9) 


Input  to  the  plant;  R  =  4sinl.884t 

Specification;  Steady  state  error  =  10%  of  average  velocity  of  the  para 
meter  drift. 

Maximum  overshoot  =  less  than  20% 

Hunting  loss  s  10%  of  optimum  parameter  value 
Expected  maximum  frequency  of  parameter  drift 
=  0.1  cps. 

Design  steps 

(1)  Gain  requirement 

Steady  state  error  =  COa 

K 

v 

where;  COa  =  average  velocity  of  parameter  drift 

K  =  static  error  coefficient 
v 

In  this  case  K  =  10 
v 

(2)  Selection  of  perturbation 
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Parameter  perturbation  signal;  sin  UJ  t 

Assume  1  v  of  parameter  adjusting  signal  corresponds  to  107o  of  optimum  para¬ 
meter  value*  Therefore,  =  1  v 

The  expected  maximum  frequency  of  parameter  drift  is  given  as  ~  O.lcps* 

Choose  CO ^  =  5  cps.  which  gives  =  1  :  50  and  gives  also  no  sharp 

overshoot  in  correction  signal  as  seen  in  simulation  study*  Thus,  =  1  sin 
31  •  At 


(3)  Establishing  a  linearized  block  diagram  of  the  system 

Referring  to  the  Fig*  4  the  complete  block  diagram  of  the  system, 
find  each  block* 

1)  G^(s)  generally,  from  3.1 


K 

Gn/S)  "  1  +  Q)  S 


The  measured  phase  characteristic  of  G  (s)  is  the  same  as  the  phase  character- 

m 

istic  of  the  plant  (Refer  to  the  section  4.1(c)(2)  ).  Therefore,  the  denom¬ 


inator  of  the  G  (s)  is  the  same  as  the  plant  transfer  function  denominator, 
m 


thus 


K 

m 

Gm(s)  =  S (1  +  0.111S) 


I  G  (s)  I  =  A  = 

I  m  I  m 


m 


(j  aO(i+o.niju;  ) 


m 


m 


A  was  measured  as  0.086  at  CO  -  0.0785  rad/sec. 


m 


is  calculated  as  0.122,  thus 


rn 


f  .  0.122 

m^S;  -  S(1+0.111S) 


2)  M. 


The  linearized  is  obtained  from  equation  (5-3) 


Gu,  =  2a  A  C, 
Ml  1  p  1 
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The  is  assumed  to  be  1.  The  attenuation  factor  of  the  perturbation 

signal  Ap  was  measured.  The  value  of  is  different  depending  on  the 

point  of  the  error  output  signal  wave  at  which  it  is  measured  as  shown  in 

Fig.  10(a). 

Maximum  A  =  0.0087 
P 

Minimum  A  =  0.00199 
P 

Average  A  =  0.0054 
P 

Take  the  maximum  A^  for  the  first  trial  design,  =  1  v,  thus 

G...  =  0.0174 
Ml 


The  value  of  A^  affects  the  overall  loop  gain  which  will  be  adjusted 


properly  at  the  final  design. 

3) 


Gl<8>  -  2 


(Vs 

1 


S“  +  +uj* 


Linearized  ^ *s>  ^ roTn  equation  (5-5), 

A 

Glequiv(s)  =  1  +  lTleq>S 


A  =  A^  =  attenuation  factor  of  G^Cs)  at  tOy  which  is  measured  as  6.6 


*  ®  ^G1(j  j  COm)  =  2. 


67' 


tan  2.67°  =  T\  S 
u  leq. 


S  =  j  CO  -  i0, 628 
m 

-x  tan  2.67' 


leq.  0.628 


=  0.0603 


Thus, 


.  .  _ 6. 6 

Gleq.  '  =  1  +  0.0603S 


«)  M2 

From  equation  (5-8),  G, 


BC. 


M2  "  2 
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B  is  the  attenuation  factor  of  reference  signal  phase  lag  filter  G^(s). 

In  this  example  G.(s)  is  not  used  because  the  G  (s)  gives  the  carrier  phase 
J  m 

shift  of  160°  which  is  closer  to  180°  phase  shift.  Therefore,  B  =  1 
Cj  s  1  v,  Thus,  G  2  =0.5 

5)  be  a  compensator  in  this  sytem  and  is  not  decided  yet. 

6)  The  integrator  G^(s)  is  chosen  to  have  following  form  because  of  diffi¬ 
culty  in  the  simulation  circuit  as  explained  in  section  4,2, 


g4(s) 


K 


4 


1  +  10S 


(4)  The  linearized  block  diagram  of  the  system  becomes  as  Fig.  35(a)  and 
which  is  reduced  to  Figure  (b) .  Excluding  GjCs)  from  the  system  and  put¬ 
ting  =  10,  the  stability  of  the  system  is  investigated.  To  obtain  = 
10,  should  be  1430.  The  open  loop  transfer  function  of  uncompensated 
system  becomes, 


Go<S>  =  S(1  +  10S)(1  +  0. Ills) <1  +  0.0603S)  ^5_9) 

The  Bode  diagram  of  the  system  is  plotted  as  Fig,  36.  The  phase  margin  is 

zero  and  the  system  is  unstable. 
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Vs) 


K, 


G,  .  (s) 
lequiv 


(a) 


0.007 iC, 


S ( 1+1CS ) ( 1+0 . 111G ) ( 1+ . 0603S ) 


(b) 


Fig*  35  The  linearized  block  diagram  of  example  1. 
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The  G^(s)  also  can  be  written  in  standard  transient  response  form  as, 


G 

o 


(s) 


_ 149 _ _ 

S(S  +  0.1)(S  +  9) (S  +  16.6) 


(5-10) 


The  root  locus  is  also  plotted  as  Fig.  37. 

The  roots  are; 

r  =  0.0336  -  jO.9874 
r2  -  0.0336  -  jO.9874 
r3  =  -9.2322 

r^  =  -16.5349  (5-11) 

Two  roots  are  located  in  the  right  half  of  the  S  plane.  Thus  the  root 
locus  plot  also  shows  the  instability  of  the  system. 

The  characteristic  equation  of  the  closed  loop  system  is, 

S4  +  25. 7S3  +  151.9S2+  14. 9S  +  K  =  0  (5-11) 

By  Routh’s  criterion*  the  critical  gain  of  the  system  is  found  as  =  87.9 
The  required  gain  of  the  system  is  149  and  exceeds  the  critical  gain. 

A  compensation  is  necessary  to  make  the  system  stable. 

(5)  Compensation 

A  lag  compensation  is  adopted,  the  compensator  is 

1  +  aTS 

g2<s>  -  TTts  <5-12) 

Compensation  1; 

To  get  5°  phase  margin  in  the  final  compensated  system,  a  cross-over  fre¬ 
quency  qj  =0.5  (phase  margin  is  7.5°  on  Bode  diagram)  is  chosen  in  Fig. 
36. 

aT  is  calculated  from  the  following  equation, 

JL  =  Jfc  (5-13) 

aT  10 
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_1  =  0^5  aT  =  20  (5-14) 

aT  10 

|Go(j^c)  |  =  12  db  on  Bode  diagram  Fig.  36 

This  must  be  attenuated,  therefore, 

20  log  a  =  -12  a  «  3-^- 


Substitute  this  value  in  equation  (5-14)  and  get  T  =  79.6.  Thus, 

q  (s\  ~  -1  +  20  S — 

G2(  }  1  +  79.6  S 


(5-15) 


This  is  the  required  compensator  to  have  5°  phase  margin. 

(6)  Investigation  of  stability  of  the  compensated  system. 

Putting  this  ^^(s)  into  the  original  system  and  you  get  the  open  loop 
transfer  function  of  the  compensated  system  as, 


p  _ 10C1+20S) _  (5-16) 

ocv }  S(1+79.6S) (1+10S) (1+0. 111S) (1+0.0603S) 

The  Bode  diagram  of  the  compensated  system  is  plotted  in  Fig.  36.  Now, 
the  system  has  5°  phase  margin  and  expected  to  be  stable. 

The  open  loop  transfer  function  of  the  compensated  system  is  rewritten 
in  standard  transient  response  form  as, 

p  ..  37 . 45  (S+0. 05) _  (5-17) 

ocV  J  “  S(S+0.0126) (S+0. 1) (S+9) (S+16.6) 

The  characteristic  equation  of  the  closed  loop  system  is 

S5+25.713S4+152.224S3+16.813S2+37.638S+1.873  =  0  (5-18) 

The  roots  are 

r  =  -0.009526  +  J0.4974 
r2  =  -0.009526  -  j0.4974 
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r*3  =  -0.05039 

ry  =  -9.0531 
4 

r5  =  -16.5905 

All  roots  are  in  the  left  half  of  the  S  plane. 

Appioximated  of  the  compensated  system  is  found  as  0.01916 

M  for  this  v/  =  1.9417 
Pt 

(7)  Simulation  of  the  compensated  system  on  analog  computer 

The  compensated  system  is  simulated  on  analog  computer.  The  actual 
response  of  the  system  is  recorded  as  Fig.  38  where  (a)  is  the  parameter 
drift  being  corrected  and  (b)  is  the  corresponding  correction  signal.  The 
response  of  the  system  is  different  depending  on  the  direction  of  drift  of 
parameter  as  explained  in  Section  4.1(b).  We  take  the  worst  case  of  response, 
that  is,  the  response  of  the  system  when  the  parameter  drifts  to  the  direc¬ 
tion  of  decreasing  system  .  From  the  recorded  result  of  the  system 

response,  following  data  are  obtained. 

Steady  state  error  =  16,57a 
M  =  approximately  2 

(8)  Comparison  of  theoretical  and  actual  data  of  response 


theoret ical  actual 


steady  state 
error 

10% 

16.5% 

M  . 

Pt 

1.9417 

2.0 
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(b)  Corresponding  correction  signal 


(b)  Example  1.1 

A  different  compensator  is  designed  for  the  example  1  and  the  compensat¬ 
ed  system  is  simulated  on  analog  computer.  The  response  of  the  actual  system 
is  also  compared  with  the  response  predicted  analytically. 

Compensation  2  ; 

(1)  To  get  10°  phase  margin  in  the  final  compensated  system,  CO =0.3 
(phase  margin  is  17°  on  Bode  diagram)  is  chosen  from  uncompensated  system 
Bode  diagram. 

Using  equation  (5-13),  aT  =  33.3 

|C0(JWc)|  =  20.5  db  (From  Bode  diagram  of  uncompensated  system) 

20  log  a  =  -20.5  a  = 

10.  O 


T  =  353 

Thus,  r  -  1  ±  33-3s 

G2l  }  1  +  353  S 


(5-19) 


This  is  the  required  compensator  to  have  10°  phase  margin  in  the  final 
compensated  system. 

(2)  Investigation  of  stability  of  the  compensated  system 

Putting  the  above  G2(s)  into  the  original  system  and  get  the  open  loop 

transfer  function  of  the  compensated  system  as, 


G  ,s) _ 10(1+33^3S) _ 

ocv  '  “  S  (1+353S)  (1+10S)  (1+0. 111S)  (1+0. 060 3S) 


(5-20) 


The  Bode  diagram  of  the  compensated  system  is  plotted  in  Fig.  39.  Now, 
the  system  has  10°  phase  margin. 

The  open  loop  transfer  function  is  rewritten  in  standard  transient  re¬ 
sponse  form  as, 


G  (s) _ 14.1  (S+0.03) _ 

oc'“  '  “  S  (S+0. 00284)  (S+0. 1)  (SJ  9)  (S+16 . 6) 


(5-21) 
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The  characteristic  equation  of  the  closed  system  is 

S5+25.70284S4+151.973S3+15.328S2+14.1423S+0,423  =  0  (5-22) 

The  roots  are 

rx  =  -0.02784  +  j0.3025 

r2  =  -0.02784  -  j0.3025 

r3  =  -0.03062 

r.  »  -9.0157 
4 

r^  =  -16.6008 

All  roots  are  in  the  left  half  of  the  S  plane. 

Approximate  of  the  compensated  system  is  0.0915 
M  for  this  J  =  1.758 

(3)  Simulation  of  the  compensated  system  on  analog  computer 

The  compensated  system  is  simulated  on  analog  computer.  The  result 
is  shown  in  Fig.  40,  where  (a)  is  the  parameter  drift  being  corrected, 

(b)  is  the  corresponding  correction  signal,  (c)  is  the  demodulated  output 
and  (d)  is  the  error  between  the  plant  output  and  model  output. 

From  the  results: 

Steady  state  error  =  11.1% 

M  =  2 

Pt 

(4)  Comparison  of  theoretical  and  actual  data  of  response 


theoret ical  actual 


steady  state 

error 

10% 

12.1% 

M 

1.758 

2.0 

Pt 

(5)  Simulation  of  compensated  system  on  digital  computer 


Because  of  difficulty  in  the  simulation  of  the  required  compensator 
on  analog  computer,  an  approximate  simulation  of  the  required  compensator 
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<d> 


1 

0 


-1 


time 


Fi£»4ft*  Example'  1.1  response  of  compensated  system 

(a)  parameter  drift  being  corrected 

(b)  Correction  signal 

(c)  Demodulated  output 

(d)  Error 


was  used  in  the  simulation  on  analog  computer  in  (3). 

The  same  system  is  simulated  on  a  digital  computer  and  the  result 
is  compared.  Refer  to  the  Appendix  I  for  the  FORTRAN  program  of  the 
simulation.  Fig.  41  is  the  graph  of  correction  signal  plotted  from  the 
numerical  data  output  of  the  computer.  If  the  Fig.  41  is  compared  with 
Fig.  40(b)  there  is  a  difference  in  signal  shape.  This  is  due  to  the  ap¬ 
proximate  simulation  of  the  compensator  in  Fig.  40. 

(c)  Example  1.2 

To  get  more  phase  margin  in  the  system  of  Example  1,  another  compensa¬ 
tor  is  designed  and  the  compensated  system  is  simulated  on  a  digital  compu¬ 
ter. 

Compensation  3: 

(1)  To  get  40°  phase  margin  in  the  final  compensated  system,  (jQ^  =  0.1 
(phase  margin  is  45°  on  Bode  diagram)  is  chosen  from  uncompensated  system 
Bode  diagram.  Calculated  in  the  same  way  as  in  example  1,  the  required  com¬ 
pensator  is  found  as 

1  +  100S 

G2(s)  =  1  +  104S 

(2)  Investigation  of  stability  of  the  compensated  system 

Putting  this  G2(s)  into  the  original  system  and  get  the  open  loop 
transfer  function  of  the  compensated  system  as, 


G 

cc 


(s) 


_ ioq+ioos) _ 

S(1+104S) (1+10S) (1+0. Ills) (1+0.0603) 


The  Bode  diagram  of  the  compensated  system  is  plotted  in  Fig.  42.  Now,  the 
system  has  40°  phase  margin. 

The  open  loop  transfer  function  is  rewritten  in  standard  transient  response 
form  as 
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G  (s)  .  -  1.49(8+0.01^ _ 

o  C  «  u. 

S(S+10  ) (S+0. 1) (S+9) (S+16. 6) 

The  characteristic  equation  of  the  closed  system  is 

S5+25.7001S^+151.90257S3+14.91519S2+1.49149S+0.0149  =  0 

The  roots  are 

rt  =  -0,04353  4-  jO, 08431 
r2  =  -0,04353  -jO, 08431 
r3  =  -0.01108 
r4  =  -8.9951 
r5  =  -16.6069 

All  roots  are  in  the  left  half  of  the  S  plane. 

From  the  pair  of  complex  roots,  of  the  system  is  found  as  0.459  and 

the  M  for  this  value  of  .f  is  1.195.  But  this  is  no  more  true  because 
these  complex  roots  are  no  more  dominant  in  the  system.  The  dominant  root 
is  r^  and  real.  Accordingly  the  expected  response  of  the  system  should  be 
an  over damped  one  and  the  time  constant  is  large. 

(3)  The  simulation  of  the  compensated  system  on  a  digital  computer 

The  compensated  system  is  simulated  on  a  digital  computer.  The  FORTRAN 
programming  for  the  simulation  on  a  digital  computer  is  shown  in  Appendix  I. 

In  calculating  for  the  simulation,  considering  an  overdamped  re¬ 
sponse.  A  =  0.00199  which  is  the  minimum  value  of  A  is  used  and  the  K, 
p  P  4 

is  found  as  6240  to  get  K  =  10.  The  result  of  the  simulation  is  shown  in 

v 

Fig.  43  which  shows  the  pole  drift  being  corrected  and  the  corresponding 
correction  signal.  The  Fig.  44  shows  the  pole  drift  being  corrected  and 
the  error.  From  the  results,  the  steady  state  error  can  not  be  found  ex¬ 
actly  because  the  signal  does  not  reach  the  steady  state  in  40  seconds  used 
in  the  computer  simulation. 
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{&)  Comparison  of  theoretical  actual  response 


It  is  predicted  theoretically  that  the  compensated  system  response  is 
overdamped  and  the  tine  constant  Is  large e  The  actual  response  shove  the 
same  pattern  of  response  as  predicted*  The  theoretical  steady  state  error 
is  lO'/L  The  aetisal  steady  state  error  c.ovild  ret  td  ™e&£M£&d  d ;/£  is  *  ^ pasted 
to  reach  close  to  the  theoretical  one* 


(d)  Example  2 


To  test  the  validity  of  the  analytical  r.eti  d  it  a-ra  extreme  case  of 


design^  the  steady  state  error  in  example  1  is  changed  to  1%  and  rigorous 
compensation  is  used0  All  other  specifications  are  the  same  as  example  1  „ 

(1)  The  required  K  «  100 

(2)  The  open  loop  transfer  function  of  rncovr.petisatea  system  is 


G 

o 


00 


_ _  100  ,,  _ 

S(1  4-  10S)(1  I  (UiiS)(i  4  0.0603~S) 


(5-23) 


The  Bode  diagram  of  uncompensated  system  Is  plotted  as  Fig,  450  The  phase 


margin  is  -29°'  and  the  system  is  unstable* 

The  G  <&)  is  also  rewritten  in  standard  transient  «^spcnse  form  as9 
o  ‘ 


G  (s) 
o 


8(5  4  0o 1) (S  *  9}{S  *  16* t) 
The  characteristic  equation  of  closed  loop  system  3s 


s4  +  25. 733  »  151. 9S*1 


14.53  +  *1 - 50  -  0 


(5*24) 


(5-251 


The  roots  arc 


r  =  0.6049  +  j  2.5461 
1 

r2  »  0.604J  -  j  2.5461 

r,  -  --11.205 
i 

t  *  -15.7003 


By  Routh v s  criterion^  the  critical  gain  is  found  as  K  ^  87,9, 

c 

The  required  gain  1490  exceeds  the  critical  gain, 

(3)  Compemsat ion 

A  lag  compensation  is  adopted.  To  get  40°  phase  margin  in  the  final 
compensated  system*  cross-over  frequency  CO^  «  0,1  {phase  margin  is  45°  on 
Bode  diagram)  is  chosen  from  uncompensated  system  Bode  diagram^  Fig ,  45 0 
Calculated  in  the  same  way  as  in  example  one  the  required  compensator  is 
found  ass 


G2(s) 


1  +  IOCS 

1  +  105s 


(5-26) 


(4)  Investigation  of  stability  of  the  compensated  system 

Putting  this  C2(s)  into  the  original  system  and  get  the  open  loop 
transfer  function  of  the  compensated  system  as9 


oc 


IS)  « 


loo(iaoos) 


(5-27) 


S(l4l0  S)<l4>10S)(l40«lllS)(U0.0603S) 

The  Bode  diagram  of  the  compensated  system  is  plotted  in  Fig,  45 ,  how* 
the  phase  margin  is  39  The  open  loop  transfer  function  is  rewritten  in 
standard  transient  response  form  as9 


r*  f+\  1,49(540,01)  (5-28) 

ocl  }  ■  irrsTcnrniJiJ^^  r, tj 

The  characteristic  equation  of  closed  system  is 

sh-25. 7000lS^»151.900257shl4.901SlS^+l.A90149St0.0149  -  C  (5-29) 

The  roots  are 

r  «  -0,04348  4-  j  0.08428 

r2  =  -0. 04348  -  JO. 08428 

r3  =  -0.01109 

r.  »  -8.935 
4 

y  5 


“  -16,6069 


All  roots  are  in  the  left  half  of  the  S  plane  and  are  the  same  values 
as  the  roots  in  the  Example  1.2.  Thus  the  same  response  as  in  the  example 
1.2  is  expected,  that  is,  the  response  of  the  compensated  system  is  over- 
damped  and  the  time  constant  is  large.  Furthermore,  a  smaller  steady  state 
error  is  expected  theoretically. 

(5)  Simulation  of  the  compensated  system  on  a  digital  computer 

The  compensated  system  is  simulated  on  a  digital  computer.  The  FORTRAN 
programming  for  the  simulation  on  a  digital  computer  is  shown  in  Appendix  1. 
In  calculating  for  the  simulation,  =  0.00199  is  used  and  the  is 
found  as  62400  to  get  =  100,  Considering  the  overdamped  response,  K, 
is  increased  by  35%  of  the  calculated  value  to  84000  in  order  to  see  the 
effect  of  an  increased  gain.  The  system  still  maintains  stability  and  the 
rise  time  is  reduced  a  little  compared  to  the  response  of  example  1.2.  The 
result  of  the  simulation  is  shown  in  Fig.  46  which  shows  the  pole  drift 
being  corrected  and  the  corresponding  correction  signal.  From  the  result, 
the  steady  state  error  can  not  be  found  because  the  signal  does  not  reach 
the  steady  state  yet. 

{6}  Comparison  of  theoretical  and  actual  response 

It  is  predicted  theoretically  that  the  compensated  system  response  is 
overdamped  and  the  time  constant  is  very  large.  The  actual  response  shows 
the  same  pattern  as  predicted  and  also  shows  that  the  35%  more  loop  gain 
than  the  designed  loop  gain  does  not  hamper  the  stability  of  the  system  at 
all,  which  tells  that  the  system  is  in  very  stable  condition  as  designed. 

The  theoretical  steady  state  error  is  1%.  The  actual  steady  state  error 
could  not  be  measured  but  the  deviation  of  the  steady  sfate  error  is  ex¬ 
pected  to  be  not  so  large. 
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6.  Conclusions. 


(a)  Analysis  of  the  system. 

The  open  loop  signal  analysis  showed  clearly  the  feasibility  of 
the  adaptive  scheme  in  this  system. 

The  signal  analysis  leads  to  the  obtaining  of  a  correct  equivalent 
signal  transfer  block  for  the  component  in  the  system. 

The  simulation  study  of  the  successive  adaptive  circuit,  from  the 
simplest  to  the  complete  circuit,  showed  the  function  of  each  component  and 
the  effect  of  component  change  in  the  circuit. 

It  is  confirmed  that  only  the  parameter  drift  signal  component  is 
acting  through  the  closed  adaptive  loop;  accordingly  only  the  parameter 
drift  signal  need  to  be  considered  in  the  closed  loop  analysis. 

The  simulation  study  showed  that  the  analytical  explanation  of  the 
system  operations  such  as  modulation,  demodulation  and  filtering,  are  in 
close  agreement  with  the  actual  operation  of  the  system, 

(b)  Design  of  the  system. 

The  obtaining  of  a  correct  error  measurement  transfer  function 

G  (s)  and  the  attenuation  factor  of  perturbation  signal  through  the  plant 

m 

(A  )  is  most  important  in  the  analytical  design  because  these  affect  the 
accuracy  of  the  analytical  design  the  greatest.  The  value  of  A^  differs 
depending  on  the  point  of  the  error  output  wave  at  which  it  is  measured. 

When  the  compensated  system  is  expected  to  be  underdamped,  high  value  of 

A  is  used  in  the  simulation  and  when  the  compensated  system  is  expected  to 

P 

be  overdamped,  low  value  of  A^  is  used  in  the  simulation  for  the  first  trial 
design. 

All  the  components  are  properly  linearized  and  an  analytical  design  is 
attempted.  The  predicted  performance  of  the  analytical  design  is  in  close 
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agreement  with  the  actual  system  performance  as  verified  in  examples. 

In  the  case  of  a  compensated  system  which  has  a  pair  of  dominant 
complex  roots  (Example  1  and  1.1)  the  steady  state  error  deviation  is 
within  107,  and  M  deviation  is  within  157*  between  the  predicted  and  actual 

pt 

value. 

In  the  case  of  a  compensated  system  in  which  the  pair  of  complex  roots 
are  not  dominant  and  a  real  root  is  dominant,  the  comparison  could  not  be 
made  numerically,  but  the  actual  time  performance  had  the  same  pattern  as 
predicted. 

This  analytical  method  of  design  simplifies  the  design  considerably 
and  is  useful  in  obtaining  a  preliminary  design  which  will  be  refined  by 
an  experimental  verification. 

It  is  further  suggested  that  this  linearized  technique  is  applied  to 
an  analysis  and  design  of  multi-loop  sinusoidal  perturbation  adaptive  con¬ 
trol  system  with  some  modifications  for  the  intercoupling  of  adaptive  loops, 
which  should  be  further  investigated. 
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APPENDIX  I 


FORTRAN  PROGRAMMING  FOR  THE  SIMULATION  OF 
A  SINUSOIDAL  PARAMETER  PERTURBATION 
ADAPTIVE  CONTROL  SYSTEM 

To  program  the  simulation  of  the  system,  the  original  block  diagram 
of  the  system  is  rearranged  to  a  proper  form  for  programming  as  shown  in 
Fig.  47  which  is  the  rearranged  block  diagram  of  the  complete  adaptive 
control  circuit,  Fig,  32.  Then,  a  simultaneous  differential  equations  are 
set  up  in  the  form  shown  at  the  end  of  the  Program  Adapt  1  (page  109  ). 

The  Program  Adapt  1  in  the  following  two  pages,  is  the  FORTRAN  program  for 
the  simulation  of  the  compensated  system  in  Example  1.1.  The  input  data 
and  the  graph  control  data  are  shown  in  proper  form  in  page  110  ,  A 
sample  numerical  output  is  shown  on  page  111. 


106 


/u,\ 

X 

1 

iV_y+r 

i 

1“ 

lo 

V: 

w 

a 

1 _ 

fr\ 

ty 

K 

i  " 

f 

rr> 

1 

X 

Ik 

.  \  \ 

C* 


107 


Fig.  47  Rearranged  block  diagram  of  a  sinusoidal  perturbation  adaptive 
control  system  for  the  simulation  on  a  digital  computer. 
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TIME 


MODEL 


PLANT 


ERROR 


CORSIG 


POLSIG 


.000 
.250 
.500 
.750 
1.000 
1.250 
1  .500 
1.750 
2.000- 
2.250- 


• OQ0OOE+OO 
.78K96E+00 
• 26836E+0 1 

•  384 1 1 E+0 1 
• 40636E+0 1 

•  34421 E+0 1 

•  207 1 3E+0 1 

•  24531 E+00 
.  16333E+01- 
• 31558E+01- 


• 00000E+03 
•65736E+03- 
• 22546E+01 ■ 

•  35108E+01  - 

•  397  44E+01  ■ 
.35493E+01 

• 23557E+01 
•64320E+03 

•  1 2369E+01 
• 28449E+01 


•OOOOOE+OO 
. 12760E+00 

•  4290  IE  +00 
•33027E+00- 

•  891  73E-0 1- 
.  10722E+00- 
.28438E+00- 
•39789E+00- 
• 39647E+00- 

•  31  088E+00- 


• OOOOOE+OO 
• 35597E-0 1 
• 80853E-01 
• 38405E+00 
• 36566E+00 
.43291 E+00 
• 36526E+00 
• 66078E+00 
. 89073E+00 
•11 579E+01 


•  OOOOOE+OO 
.  60356E+01 
•52317E+01 
• 36684E+01 

•  44477E+01 
•55669E+01 
•47988E+01 

•  33968E+01 

•  39046E  +  01 
•484 1 4E+01 


2.500- 
2.750- 
3.000- 
3.250- 
3.  500 

3.750 
4.000 
4.250 
4.500 

4.750 


•  399 1 OE+O 1 
• 39577E+0 1 
. 30625E+0 1 
. 14998E+01 
• 38917E+00 
• 21934E+0 
• 35204E+0 
•40801E+0 
.37512E+0 
.26061  E+0 


•381 31E+01 
.39432E+01 
••32066E+01 
■•17681  E+01 
•51 772E-01 

•  1 874 1 E+01 
•32972E+01 

•  3981  5E+01 
•37897E+01 
•27777E+01 


. 17783E+00- 
.  14467E-01- 
14408E+00- 
•26829E+00- 
’•  33740E+00- 
31935E+00- 
’.22322E+00- 
-.98508E-01- 
• 38573E-01- 
. 17165E+00- 


. 13639E+01 
.  12966E+01 
.1 3352E+01 
•1351 5E+01 
.  14774E  +  01 
.171 10E+01 
•  1 781 2E  +  01 
.1 8987E+01 
. 1 8800E+0 1 
.  18804E+0 1 


.  381 35E+01 
•27664E+01 
•3442 1E+01 
•  46470E+01 
•3713 1E+01 
.23575E+01 
•29782E+01 
• 40989E+01 
• 33238E+01 
•21 939E+01 


5.000  . 
5.250-. 

5.500- . 

5.750- . 
6.000-. 
6* 250—. 

6.500- . 

6.750- . 
7.000  . 
7.250  . 


89316E+00 
10144E+01- 
27007E+01- 
37993E+01- 
40708E+01- 
34550E+0  1- 
20869E  +  0 1- 
26484E+00- 
1 61  53  E+0  1 
31 434E+0 1 


. 1 1495E+01 
■  73698E+03 
•24559E+01 

•  36401 E+01 
.40252E+01 

. 35242E+01  ■ 
•22648E+01 

•  5 1 085E+03- 

•  137  62E+01  ■ 

•  2  9605E+01  • 


•  25633E+00- 

•  27742E+00- 

•  24480E+00- 
.1  5925E+00- 
. 45597E-01  ■ 
.6921  OE-0 1- 
>•  1  7786E+00- 
.24601E+00- 
>•  2391 OE+OO- 
•.  18291E+00- 


-.201  19E+01 
•  2  1 224E+0 1 
■.231  52E+01 
■•23608E+01 
■•23360E+01 
■•23203E+01 
••  22784E+01 
•.24013E+01 
■•25276E+01 
.26517E+01 


. 27297E+01 
• 38739E+01 

•  290 1 7E  +01 

•  1 7 1 94  E+01 
• 23879E+01 
• 36745E+01 
•2951 7E+01 

•  1  685 1 E  +01 

•  2 1 787E  +  01 
•3341 3E+01 


7.500 

7.750 

8.000 

8.250 

8.500- 

8.750- 
9.000- 
9.250- 

9.500- 

9.750- 


• 39867E+0 1 
. 39625E+0 1 
. 30754E+01 
•151 80E+0 1 

•  36972  E+00" 
.21769E+01- 

•  35 1 04E+0  1- 

•  40788E+0  1- 
.  37589E+01- 

•  262 1 OE+0 1- 


•  38835E+01  ■ 
.39601  E+01  ■ 
.3  17  17  E+01 

•  1 6909E+01 
*•  1  5033E+0D 
•.1  9676E+01 
’•33663E+01 
-.401  72E+01 
’•  37887E+01  ■ 
■  .27431  E+01  ■ 


•.  10322E+00- 
-.2361  9E-02- 

•  96302E-01- 

•  1 7297E+00- 

•  2 1 939E+00- 

•  20928E+00- 
.  14413E+00- 

•  61  639E-0 1* 

•  29821 E-Ol  ■ 
.  1  2207E+00- 


■.272896+01 
•  • 267  58E+0 1 
26633E  +  01 
■•26638E+01 
■.27329E+01 
•28556E+01 
28980E  +  01 
■  • 29286E+01 
•.28911  E+01 
.28671 E+01 


.  251 42E+01 
.  14171  E+01 
• 20254E+01 
• 33272E+01 
. 25233E+01 
. 1 2438  E+01 
. 17734E+01 
• 30600E+01 
.237816+01 
. 12391E+01 


10.000- 

10.250 

10.500 

10.750 
1  1  .000 

11.250 

11.500 

1 1.750 
12.000- 
12.250- 


■  91 222E+00- 
• 99548E+00 
• 26860E+0 1 
•37920E+0 1 
• 40725E+01 
• 34654E+0 1 
.210376+01 
•28433E+00 
. 159736+01 
.31309E+01- 


-.  10906E+01 
•80536E+00 
•251 65E+01 

•  36806E+01 
.404  18E+01 

•  35 1 62E+01 
•22350E+01 
.46815E+00 

-.142026+01 

-.299706+01 


-.  17834E+00- 
-.  1901 1E+00- 
-.  16951E+00- 
-.11  141E+00- 
-.30716E-01- 

•  50821  E-Ol- 

•  1 3 1 30E+00- 
.183826+00- 

•  1 771  5E+00- 

•  13390E+00- 


.292986+01 
• 29963E+01 
.310746+01 
•31 278E+01 
•3095 1 E+0 1 
•30610E+01 
.302246+01 
•30866E+01 
.317136+01 
.324186+01 


.  17242E+01 
•29897E+01 
• 21748E+01 
.985396+00 
.  15418E+01 
•29220E+01 
• 22727E+01 
.103376+01 
.  14486E+01 
.27381 E+01 
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